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ABSTRACT 


This  report  summarizes  a  body  of  research  dealing  with  the  nonlinear  control  theory 
of  complex  mechanical  systems.  The  principal  focus  is  on  so-called  super- articulated  (or 
underactuated)  mechanical  systems — systems  having  a  mixture  of  directly  controlled  as  well 
as  uncontrolled  degrees  of  freedom.  The  twofold  goal  of  the  research  is:  (1)  how  to  control 
the  degrees  of  freedom  which  are  directly  actuated  without  eliciting  undesired  behavior  in  the 
unactuated  degrees  of  freedom,  and  (2)  how  to  prescribe  motions  of  the  directly  actuated 
degrees  of  freedom  to  achieve  motion  objectives  for  the  degrees  of  freedom  which  are  not 
directly  controlled. 

The  report  contains  a  comprehensive  list  of  research  published  under  this  grant  as  well  as 
some  predecessor  grants.  Two  papers  in  particular — “The  Geometry  of  Controlled  Mechanical 
Systems”  and  "Open-loop  Oscillatory  Stabilization  of  an  n-Pendulum”  comprise  the  present 
report. 


The  Geometry  of  Controlled  Mechanical  Systems 

J.  Baillieul* 

April  13,  1998 


Abstract 

The  control  of  mechanical  systems  has  become  a  principal  application  focus  of  nonlinear 
control  theory.  This  development  was  to  a  very  large  extent  foreseen  in  the  early  work  of 
R,W,  Brockett,  and  a  number  of  his  papers  from  the  mid  1970’s  have  been  highly  influential 
in  establishing  the  links  between  nonlinear  control  theory  and  geometric  mechanics.  Taking 
some  of  Brockett’s  early  work  as  a  starting  point,  we  study  intrinsically  second  order  nonlinear 
control  systems.  The  theory  is  developed  using  the  language  of  differential  geometry  and  affine 
connections.  Velocity  and  acceleration  controlled  Lagrangian  systems  provide  a  rich  class  of 
examples  of  second  order  control  systems.  These  cirise  in  modeling  the  dynamics  of  super- 
articulated  mechanical  systems,  in  which  only  some  of  the  configuration  variables  (or  degrees 
of  freedom)  are  directly  controlled,  with  the  remaining  variables  evolving  under  the  dynamic 
influence  of  the  actuated  degrees  of  freedom.  Our  goal  is  to  develop  a  nonlinear  control  theory 
which  characterizes  the  way  in  which  the  unactuated  degrees  of  freedom  in  these  systems  are 
influenced  by  the  degrees  of  freedom  which  can  be  controlled  directly.  Central  to  the  theory  is  a 
natural  Riemannian  structure  and  certain  curvature-like  quantities  which  provide  a  measure  of 
the  lack  of  integrability  in  the  mapping  of  input  trajectories  to  configuration  space  trajectories. 
The  report  chapter  concludes  with  results  connecting  the  geometry  of  second  order  systems  with 
the  theory  of  averaged  Lagrangian  and  Hamiltonian  systems. 


1  Introduction 

The  early  years  (=for  the  most  part  the  decade  of  the  1970’s)  in  the  development  of  nonlinear  geo¬ 
metric  control  theory  witnessed  attempts  to  extend  all  of  what  was  known  about  linear  systems  into 
the  realm  of  controlled  dynamical  systems  which  could  be  described  by  a  finite  number  of  ordinary 
differential  equations  with  inputs.  A  number  of  researchers  joined  the  quest,  and  a  sampling  of  sem¬ 
inal  papers  includes  Brockett,  [1972,  1973],  Haynes  and  Hermes,  [1970],  Hermes  and  Haynes,  [1963], 
Krener,  [1974],  and  Sussmann  and  Jurdjevic,  [1972].  The  simple  observation  that  differentiable 
manifolds  were  the  natural  state  space  structures  for  the  study  of  controlled  nonlinear  differential 
equations  opened  a  field  of  research  that  would  be  rich  and  active  for  many  years. 

While  the  field  developed  initially  with  the  goal  of  extending  linear  state  space  control  theory 
into  the  nonlinear  domain,  there  was  the  early  discovery  of  features  of  nonlinear  systems  for  which 
no  linear  analogues  exist.  Specifically,  the  geometry  of  noncommuting  vectorfields  is  a  theme  which 
was  central  in  the  earliest  expositions  of  geometric  nonlinear  control  theory,  and  this  theme  has 
remained  central  to  many  major  developments  over  the  past  quarter  century.  In  addition  to  laying 
down  significant  parts  of  the  theoretical  foundations  of  modern  geometric  control  theory,  Roger 
Brockett  also  did  seminal  work  on  important  applications  of  the  theory  to  the  study  of  controlled 

*The  author  is  very  grateful  to  Dr.  Marc  Jacobs  of  the  Air  Force  Office  of  Scientific  Research  for  many  years  of 
support  of  his  research  on  the  nonlinear  control  theory  of  mechanical  systems.  Currently  this  support  is  extended 
through  grant  number  F49620-96- 1-0059.  Support  from  the  Army  Research  Office  under  the  ODDR&E  MURI97 
Program  Grant  No.  DAAG55-97-1-0114  to  the  Center  for  Dynamics  and  Control  of  Smart  Structures  (through  Harvard 
University)  is  also  gratefully  acknowledged.  Finally,  it  is  a  pleasure  to  thank  P.S.  Krishnaprasad  for  his  insights  on 
the  matters  discussed  herein. 
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mechanical  systems.  The  1976  NASA-Ames  lectures  (Brockett,  [1976])  began  a  school  of  research 
that  has  now  grown  too  large  to  survey  in  the  present  brief  article.  It  is  probably  not  coincidental 
that  during  the  twenty-five  year  development  of  geometric  nonlinear  control  theory,  the  field  of 
analytical  mechanics  has  itself  been  developed  along  increasingly  geometric  lines.  (See,  for  instance 
Marsden  and  Ratiu,  [1994].) 

The  starting  point  for  this  article  is  the  global  definition  of  a  smooth  nonlinear  control  system 
which  appeared  for  the  first  time  in  Brockett,  [1976],  In  the  next  section,  we  repeat  this  definition  and 
discuss  its  extension  and  specialization  to  the  case  of  intrinsically  second  order  systems.  While  such 
systems  can  be  defined  and  studied  in  a  purely  formal  way,  we  shall  show  that  they  arise  naturally 
as  Lagrangian  models  of  super- articulated  (also  referred  to  as  under- actuated)  mechanical  systems. 
In  these  systems,  some  of  the  configuration  variables  (or  degrees  of  freedom)  are  directly  controlled 
while  the  remaining  variables  evolve  under  the  dynamic  influence  of  the  controlled  degrees  of  freedom. 
We  shall  study  such  mechanical  systems  taking  the  viewpoint  that  the  directly  actuated  degrees  of 
freedom  are  inputs  or  controls  which  influence  the  remaining  configuration  variables  through  the 
given  Lagrangian  dynamics.  Of  course,  in  Lagrangian  models,  the  velocities  and  accelerations  of  the 
directly  actuated  degrees  of  freedom  also  play  a  role  in  influencing  the  dynamics  of  the  unactuated 
degrees  of  freedom,  and  thus  it  is  natural  to  regard  the  inputs  to  these  systems  as  triples  comprised 
of  the  actuated  configuration  variables  together  with  their  associated  velocities  and  accelerations. 

In  Section  3,  we  shall  associate  certain  curvature-like  quantities  to  our  super-articulated  La¬ 
grangian  control  systems.  Both  configuration  flatness  and  input  flatness  will  be  defined,  and  it  will 
be  shown  that  input  flatness  and  configuration  flatness  together  are  necessary  conditions  for  there  to 
exist  a  change  of  coordinates  such  that  in  the  new  coordinates  (i)  the  inertia  tensor  is  diagonal,  and 
(a)  no  input  accelerations  enter  the  equations  of  motion.  (In  this  case,  we  may  say  the  systems  are 
velocity  controlled.)  It  will  be  noted  that  the  flatness  discussed  in  Section  3  is  an  intrinsically  second 
order  concept — related  to  ideas  from  classical  Riemannian  geometry.  Thus  we  must  distinguish  this 
from  differential  flatness  which  has  been  introduced  by  Fliess  (e.g.  [15],  and  see  also  [23]). 

In  Section  4,  we  consider  the  same  Lagrangian  models  of  super-articulated  mechanical  systems 
under  the  influence  of  oscillatory  inputs.  For  this  class  of  systems,  we  wish  to  answer  the  question  of 
when  high  frequency  forcing  of  certain  degrees  of  freedom  will  cause  the  total  assemblage  to  respond 
with  a  prescribed  stable  behavior.  A  growing  body  of  research  on  this  question  has  been  aimed 
at  understanding  the  extent  to  which  this  question  can  be  answered  by  an  appropriate  theory  of 
averaging.  Section  4  presents  recent  results  on  averaging  Hamiltonian  and  Lagrangian  systems,  and 
in  particular,  we  describe  the  role  of  the  averaged  potential  in  characterizing  the  qualitative  dynamics 
of  periodically  forced  mechanical  systems.  Configuration  flatness  together  with  input  flatness  again 
play  a  role  in  characterizing  necessary  conditions  for  the  existence  of  a  coordinate  system  in  terms  of 
which  the  averaged  potential  is  obtained  directly  by  simple  averaging  of  the  potential  energy  terms 
in  the  Lagrangian. 

2  Second  Order  Generalized  Control  Systems 

A  local  description  of  a  nonlinear  control  system  typically  involves  a  differential  equation  of  the  form 

q  =  fiq,x),  (1) 

where  f  :  M  xU  TM,  with  M  being  an  (n-dimensional)  state  manifold  and  U  an  (m-dimensional) 
input  manifold.  Brockett  [1976]  pointed  out  that  such  models  do  not  provide  global  descriptions  of 
systems  where  the  available  input  is  state  dependent.  In  the  natural  global  extension  of  (1),  we  view 
/  as  a  bundle  map. 

Definition  1  A  smooth  nonlinear  control  system  is  a  quadruple  (B,M,7r, /)  such  that 

(i)  (jB,  M,  tt)  is  a  fiber  bundle  with  total  space  B,  base  space  M,  and  canonical  projection  tt  :  B  M, 
and 
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(ii)  f  :  B  TM  is  a  bundle  morphism  such  that  for  each  q  £  M  and  x  £  Uq  =  t:  ^{q), 
f{q,x)  £  TqM. 

For  a  complete  introduction  to  the  theory  of  fiber  bundles,  the  reader  is  referred  to  Husemoller 
[1994]  or  the  classical  treatise  Steenrod  [1951].  Recall  that  for  each  q  £  there  is  a  neighborhood 
V  of  q  and  a  diffeomorphism  (j)  mapping  7r”^(y)  onto  V  x  7r^^(g).  Thus,  as  Brockett  points  out 
(Brockett,  [1976],  p.  16),  by  restricting  our  attention  to  such  neighborhoods,  it  is  always  possible  to 
find  a  local  representation  of  a  smooth  nonlinear  control  system  which  is  of  the  form  (1). 

To  study  controlled  mechanical  systems,  it  will  be  useful  to  extend  models  of  the  form  (1)  to 
include  systems  which  are  second  order  (in  the  configuration  variables  q)  and  in  which  the  controlling 
effects  of  the  input  variables  are  primarily  due  to  their  accelerations,  ^(*).  The  descriptions  of  motions 
in  these  systems  will  involve  the  notion  of  parallel  displacements  of  vectorfields  in  TM  along  curves 
in  B.  Given  a  vectorfield  X  on  B  and  any  curve  r(t)  £  B,  the  lifting  of  r  to  a  curve  X{r{t))  £TB 
is  projected  by  the  tangent  mapping  tt*  onto  a  curve  7r*(X(r(t)))  £  TM.  This  association  leads 
to  the  desired  notion  of  covariant  differentiation  of  vectorfields  in  TM  along  curves  in  B,  and  we 
shall  prescribe  this  in  terms  of  local  coordinates  as  follows.  Let  (gi , . . . ,  2:1 , . . . ,  Xm)  denote  local 

coordinates  defined  in  a  neighborhood  of  a  point  p  £  B.  Let 


II 

.,n), 

II 

Ii 

..,m) 

be  the  associated  vectorfields  in  TB.  We  may  also  view  (gi, . . .  ,gn)  as  defining  local  coordinates  in 
a  neighborhood  of  7r(p)  £  M.  Then  we  similarly  let 

=  . ") 

be  the  associated  vectorfields  in  TM.  In  terms  of  this  local  coordinate  description,  for  each  Xi  and 
each  (g,  x)  £  B,  the  tangent  mapping  tt*  associates  a  vector  Y^*  £  TqM  by  means  of  the  formula 

F;(g,a:)=7r*Xi(g,a:). 


In  general,  y/(g,x)  7^  Yiiq),  but  we  may  write 

n 

=  Y (*  =  1. •••.«)• 

i=i 


By  restricting  to  a  smaller  neighborhood  if  necessary,  there  is  no  loss  of  generality  in  assuming  that 
the  nxn  matrix  A{x),  whose  ij-th.  entry  is  aij(x),  is  nonsingular  for  each  x.  For  any  point  q  £  M, 
^  each  vectorfield  Xi  {i  =  l,...,n)  and  Xj  {j  =  defines  an  integral  curve 

passing  through  7r"’^(g).  For  each  such  integral  curve,  expressed  in  our  coordinate  neighborhood  as 
r{t)  =  {q{t)^x{t))^  there  is  a  corresponding  set  of  curves  Y*{q{t),x{t))  (i  =  1 . . .  ,n)  in  the  tangent 
bundle  TM.  Tangents  to  these  curves  define  a  covariant  derivative  operator  V  :TB  xTM  -¥TM 
which  is  prescribed  in  terms  of  the  given  local  coordinates  by  means  of  a  total  of  (n+m)n^  functions: 
{r^}  {i,j,  fc  =  1, . . . ,  n)  and  (i  =  1, . . .  ,m\  j,k  =  1,. n),  and  the  formulas 


^x,y;  = 
= 


En  p/.  y* 

k=l  ^  k  > 

E”  pA;  V* 


z  =  1, . . .  ,m; 


(2) 

(3) 


Since  the  vectorfields  Xi  (and  Xi)  lie  in  the  coordinate  directions,  it  is  clear  what  is  meant  if  we 
write 


(4) 


k=l 


k=l 
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The  covariant  derivative  operator  defines  a  generalized  affine  connection  and  allows  us  to  study 
displacement  of  vectorfields  in  TM  along  curves  {q{t)jx{t))  G  B.  We  refer  to  Kobayashi  and  No- 
mizu,  [19],  for  a  classical  treatment  of  affine  connections  and  parallel  displacement.  In  the  present 
development,  we  shall  view  points  q  in  the  base  manifold  as  states  of  a  nonlinear  control  system  (and 
eventually  as  generalized  coordinates  of  a  mechanical  control  system)  and  the  points  x  in  the  fiber 
as  control  inputs.  We  shall  be  interested  in  lifting  curves  in  the  total  space  B  to  curves  in  TB  and 
projecting  these  liftings  onto  TM.  It  is  useful  to  keep  in  mind  the  following  commutative  diagram 

TB  TM 

TT 

B  - ^  M 


where  the  vertical  arrows  are  the  canonical  tangent  bundle  projections.  As  above,  given  a  curve 
{q{t),x{t))  in  B,  we  lift  it  to  a  family  of  curves  Xi{q{t),x{t))  in  TB  {1  <  i  <  n)  and  define  corre¬ 
sponding  projections  Y*{q{t),x{t))  =  TT*Xi{q{t),x{t))  in  TM.  If  {q{t),x(t))  {0  <t  <  to)  remains  in 
our  coordinate  neighborhood  in  B,  we  have  a  valid  local  representation  of  the  corresponding  tangent 
vector 

n  m 

ai(t)Xi(g(t),x(t))  +  Mt)Xi(g(t),x(t)) 

i=l  t=l 

for  each  t  in  0  <  t  <  to.  (This  is  a  curve  in  TB.)  The  vectorfields  y]*(qfX)  =  7r*Xi(gf,x),  appearing 
in  the  formulas  (2)-(3)  thus  depend  on  points  in  the  total  space  B.  In  terms  of  the  vectorfields  Y*, 
the  image  curve  in  TM  may  be  expressed  as 

1=1 

This  projection  defines  a  vectorfield  (in  TM)  along  the  curve  T{t)  =  {q{t),x{t))  in  B.  Differentiating 
Y  with  respect  to  t  and  using  our  covariant  derivative  operator  to  project  the  result  back  onto  TM, 
we  obtain  the  covariant  derivative  in  the  direction  f,  for  each  t  in  0  <  t  <  to: 

n  /  n  m  \ 

i=i  ''  j=i  j=i  ' 

(Cf.  Kobayashi  and  Nomizu,  [19]  p.  114.)  Using  the  relationships  (4),  this  expression  may  be  rendered 

n  y  n  n  m  \ 

E  +  E  +  EEf.Wi)n*-  (6) 

k=l  '  *>j=l  ^ 

Using  the  notion  of  parallel  displacement^  we  can  now  state  conditions  under  which  the  curves  of  the 
form  (5)  are  trajectories  of  a  second  order  generalized  control  system.  The  conditions  involve  the 
vanishing  of  covariant  derivatives,  and  thus  these  curves  are  natural  generalizations  to  the  input- 
dependent  case  of  geodesic  flows. 

Definition  2  We  say  that  a  curve 


i=l 
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is  parallel  along  T{t)  =  if  the  covariant  derivative  in  the  direction  f,  V-f-V  is  zero.  Equiv¬ 

alently,  in  terms  of  local  coordinates,  V{t)  is  parallel  along  r(t)  if 


m  V 


=  0.D 


An  even  more  explicit  rendering  of  the  condition  for  Y  in  (5)  to  be  parallel  along  r  is  that  a(*) 
is  the  unique  solution  to  the  system  of  dififerential  equations 


n  n  m 

+  E  +  E  E =  0  (A:  =  1, . . . , n). 

ij=l  i=l  j=l 

Definition  3  We  say  that  (5)  defines  a  smooth  generalized  second  order  control  system  if  (i)  ai{t)  = 
+  Yl]Li  7ii(9)^j(t)  where  the  'jij  are  mn  smooth  functions  of  q,  and  (ii)  the  vector  field  (5)  is 
parallel  along  the  curve  {q{t),x{t)).D 


In  particular,  if  (5)  is  a  trajectory  of  smooth  generalized  second  order  control  system,  the  coor¬ 
dinate  functions  qi{t)  and  Xi{t)  satisfy 


m  n  n  m 

%  +  E  'Tij  (9)^i  E  =  0,  for  i  =  1, . . . ,  n,  (7) 

jz=l  ij=l  t=l  i=l 


where  f  =  fj^  +  This  is  a  direct  extension  of  the  geometric  characterization  of  second-order 
differential  equations  found  in  Abraham  and  Marsden,  1988,  [1].  We  employ  the  word  generalized  to 
indicate  that  derivatives  of  the  input  curves  x{')  influence  the  dynamics  in  a  direct  and  significant 
way.  (See  Fliess,  [1990],  for  a  broader  discussion  of  generalized  control  systems.)  In  general,  all 
three  components  of  the  triple  (x,  x,x)  enter  the  equation  (7),  but  in  Sections  4  and  5  we  shall 
be  principally  concerned  with  high-frequency  periodic  inputs,  and  there  the  most  significant  terms 
will  be  those  involving  x.  (Suppose  x(‘)  is  a  periodic  vector  valued  function  whose  fundamental 
frequency  is  lj  and  whose  £oo  norm  is  0{1).  Then  the  norm  of  x(-)  is  0{uj)  and  the  norm  of  x(-)  is 
0{u;^).)  When  equation  (7)  depends  on  x  that  we  shall  call  these  systems  acceleration  controlled. 
In  the  next  section,  we  shall  study  systems  in  which  all  jij{q)  =  0.  We  shall  also  show  that  the 
vanishing  of  curvature-like  quantities  is  a  necessary  condition  for  there  to  be  a  transformation  to 
coordinates  in  terms  of  which  x  and  x  (but  not  x)  enter  the  equations  of  motion. 

Generalized  Lagrangian  Control  Systems  Generalized  second  order  control  systems  arise 
in  the  study  of  super- articulated  mechanical  systems.  (See  Seto  and  Baillieul  [1994]  for  the  basic 
theory.  These  have  also  been  called  under- actuated  systems  in  the  literature.)  The  general  framework 
assumes  there  is  given  a  Lagrangian  L{y,y)  defined  on  the  tangent  bundle  TQ  of  the  configuration 
space  of  a  mechanical  system.  Suppose  that  the  generalized  coordinates  can  be  partitioned  as  y  = 
(r,  q)  and  that  exogenous  generalized  forces  (control  inputs)  can  be  applied  to  only  the  coordinates  r, 
while  the  coordinate  variables  comprising  q  evolve  freely,  subject  only  to  dynamic  interactions  with 
the  r-variables.  In  this  partitioning,  r  is  an  m-tuple  and  q  is  an  n-tuple  of  generalized  coordinate 
variables.  The  equations  of  motion  for  the  system  take  the  form 


dt  dr  dr 


(8) 


dt  dq  dq  ’ 


(9) 


where  u  is  an  m-vector  of  controls,  and  we  suppose  the  mapping  q 
if  the  Lagrangian  takes  the  form  “kinetic  minus  potential  energy,”  i.e. 


the  inertia  matrix  is  partitioned  conformably  with  (r,g),  M  = 


^  is  invertible.  That  is, 
L  =  ^y'^My  -  V{y),  and 

then  M  =  -M(r, g) 
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is  an  n  X  n  invertible  matrix.  It  is  also  assumed  that  M  =  A/’(r,  q)  is  an  invertible  matrix,  and 
this  ensures  a  one-one  correspondence  between  control  trajectories  u(*)  and  trajectories  of  of  the 
r- variables.  For  the  purposes  of  the  model,  we  assume  that  the  components  Ui(*)  are  each  piecewise 
analytic  functions  on  [0,oo),  although  more  general  inputs  (e.g.  impulse  trains)  are  also  of  interest 
and  amenable  to  study. 

A  simple  reduction  process  eliminating  the  explicit  dependence  on  the  input  u(-)  leads  to  a  system 
of  equations  of  the  desired  form  in  which  the  r-variables  and  their  derivatives  play  the  role  of  inputs 
with  the  triple  {x^v,a)  =  (r,r,r)  controlling  the  state  variables  (g,g)  according  to  equation  (9). 
In  this  case  the  dynamical  relationship  (9)  may  also  be  obtained  by  applying  the  Euler-Lagrange 
operator  ^  ^  ^  to  tho  reduced  Lagarangian 

C{q,q]X,v)  =  ^q^ M{q,x)q  +  v'^Aiq,x)q  -  V{q\x,v),  (10) 

where  V(g;a;,u)  =  V{q,x)  —  ^v'^J^{qyx)v.  Writing  out  the  equations  of  motion  explicitly,  we  have 

n  m  n  n  m 

^  ^  ^  ^  "b  ^  ^  ^ijkQiQj  “b  ^  ^  ^  ^  ^Cjk'^CQj  —  1, . .  . ,  n),  (11) 

i=l  ^=1  i,j=l  j=l 

where 

_l/dmki  orrikj  dmij\ 

'^^^2\dqj  ^  dqi  dqk)' 
p  _  drrikj  datk  daij 

dxt  %  dq,^ 

and  aij  and  mij  are  the  ij-th  entries  in  the  mxn  and  nxn  matrices  A{q,  x)  and  M(q,  x)  respectively. 
F{t)  is  a  vector  of  generalized  forces  Fi{t)  =  These  may  be  thought  of  as 

coming  from  a  velocity-dependent  potential  if 

d^ati  _  d^aij 
dqjdxk  dqidxk 

for  all  =  1, . . .  ,m  and  =  1 . . .  ,n.  In  classical  mechanics,  the  quantities  defined  in  this 
way  in  terms  of  the  inertia  tensor  M  are  called  Christoffel  symbols  of  the  first  kind.  To  be  consistent 
with  this  nomenclature,  we  call  the  tijk  input  symbols  of  the  first  kind. 

To  define  corresponding  symbols  of  the  second  kind,  let  denote  the  ij—th.  element  of  . 
Multiplying  both  sides  of  (11)  by  and  summing  over  the  index  values  fc  =  1, . . .  ,n,  we  obtain 

m  n  n  m 

Q(T  "b  ^  ^  H-  ^  ^  ^ijQiQj  "b  V  ^  ^  ^  “  ■^(t(0>  ((J  =  1,  .  .  .  ,  Tl)^  (12) 

£=1  ij=l  j=l  ^=1 

where 

n 

Jat  =  X) 

fc=l 

k=l 

n 

%  =  E  and 

k=\ 

n 

F„{t)  =  Y. 

*!=1 
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f  The  quantities  F^-  are  called  Christoffel  symbols  of  the  second 

kind^  and  the  FJj’s  will  be  called  input  symbols  of  the  second  kind.  Modulo  a  minor  change  of 
notation  {xj  =  vj),  when  the  generalize  potential  forces  Fk{t)  =  0,  equations  (7)  and  (12)  are  the 
same.  To  be  consistent  with  the  terminology  introduced  above  for  abstract  second  order  systems, 
we  shall  call  systems  of  the  form  (12)  generalized  Lagrangian  control  systems.  Whenever  at  least 
one  jij  ^  0,  i  =  l,...,n;  j  =  l,...,m,  these  could  also  be  referred  as  acceleration- controlled 
Lagrangian  systems.  This  usage  is  quite  informal  at  this  point,  and  it  will  be  shown  below  that 
under  certain  geometric  conditions  there  is  a  change  of  coordinates  which  eliminates  the  explicit 
dependence  of  the  dynamics  on  any  of  the  vi  terms.  Indeed,  one  of  the  main  results  of  the  chapter 
is  to  provide  an  invariant  characterization  of  generalized  Lagrangian  control  systems  which  are 
intrinsically  acceleration  controlled. 

Remark  1  Connections  with  classical  Lagrangian  reduction.  Marsden  [1992]  describes  Lagrangian 
reduction  for  mechanical  systems  with  symmetries  which  are  characterized  in  terms  of  a  group 
action  that  leaves  the  (unreduced)  Lagrangian  invariant.  For  Abelian  group  actions,  this  reduction 
was  known  to  Routh  in  the  nineteenth  century,  but  it  has  been  extended  to  the  non-Abelian  case 
by  Marsden  and  various  others  cited  in  Marsden  [1992].  The  idea  is  that  when  a  Lagrangian  is 
constant  on  orbits  of  the  tangent  lift  of  a  group  action  by  a  symmetry  group  G  acting  freely  on 
the  configuration  manifold  Q,  these  orbits  can  be  “factored  out”  to  yield  a  reduced  configuration 
space  Q/G  on  which  there  is  defined  a  reduced  Lagrangian  system.  Symmetries  are  the  key  to 
this  reduction.  In  the  above  reduction  of  Lagrangian  systems  with  inputs,  we  have  assumed  no 
particular  structure  other  than  a  natural  partitioning  of  configuration  space  variable  into  two  sets: 
those  variables  that  are  directly  controlled  and  those  variables  whose  motions  arise  solely  from 
dynamic  interactions.  Without  further  assumptions,  there  is  little  that  can  be  said  in  general  about 
the  dynamics  of  our  reduced  Lagrangian  systems  (12).  In  the  next  section,  we  shall  show  that 
curvature-like  quantities,  defined  by  the  Christoffel  symbols  and  input  symbols,  and  related  notions 
of  flatness  play  a  role  in  characterizing  key  features  of  systems  of  the  form  (12).  When  these 
quantities  vanish,  there  is  a  choice  of  coordinates  with  respect  to  which  the  inputs  and  configuration 
variables  are  highly  decoupled.  The  “vanishing”  of  certain  additional  curvature-like  quantities  is  a 
necessary  condition  for  there  to  be  a  choice  of  coordinates  such  that  the  dynamic  effects  of  input 
triples  {x,x,x)  are  completely  determined  by  the  configuration-velocity  pair  {x,x)  alone.  These 
terms  will  enter  the  Lagrangian  model  through  velocity  dependent  potential  terms. 

Remark  2  Invariance  properties  of  the  reduced  Lagrangian.  An  important  feature  of  Lagrangian 
mechanics  is  the  invariance  of  Lagrange’s  equations  with  respect  to  coordinate  changes.  Reduced 
Lagrangians  of  the  form  (10)  enjoy  essentially  the  same  invariance  with  respect  to  coordinate  trans¬ 
formations  of  the  form  Q  =  F{q)  which  are  independent  of  the  input  x.  Letting  the  inverse  trans¬ 
formation  be  written  q  =  G(Q),  the  reduced  Lagrangian  in  Q-coordinates  is  expressed  as 

C{Q,Q]x,v)  =  M{Q,x)Q  +  v'^ A{Q,x)Q  -V{Q\x^v),  (13) 

where 

M{Q,x)  =  ^  M{G{Q),x)-^ 

A(Q,x)  =  A{G{Q),x)^,  and 

V(Q;x,ii)  =  V(G(Q);x,t;). 

Invariance  with  respect  to  such  a  coordinate  transformation  means  the  equations  of  motion  for  the 
system  in  Q-coordinates  are  given  by 

dt  dQ  dQ~  ' 
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In  the  following  sections,  it  will  be  of  interest  to  allow  coordinate  transformations  of  the  the  config¬ 
uration  variables  which  are  input  dependent:  Q  =  F(g,x).  In  this  case  it  is  possible  to  propose  an 
extended  notion  of  invariance  for  the  equations  of  motion.  Specifically,  suppose  that  for  some  (go,  a:o) 
there  is  a  neighborhood  U  xV  CiC  'x  H""  containing  (go,a;o)  and  a  mapping  F  :U  xV 
such  that  for  each  x  E  F,  Q  =  F(g,x)  is  a  diffeomorphism  on  U.  For  each  a:,  write  the  inverse  of 
this  configuration  space  diffeomorphism  as  g  =  G{Q,x),  Then  in  terms  of  the  Q-coordinates,  the 
reduced  Lagrangian  again  has  the  form  (13)  where 

MiQ,x)  =  ^  M(G(Q,x),x)  — 

A{Q,x)  =  MiG{Q,x),x)  ■^+A{G{Q,x),x)  and 

V{Q;x,v)  =  -x'^(^A{G{Q,x),x)  M{G{Q,x),x)  |f)x 

+V((?((5,x);x,t>). 


(14) 

(15) 

(16) 


Again,  we  find  that  our  reduced  Lagrangian  formulation  is  invariant  under  this  class  of  input 
dependent  coordinate  transformations.  This  is  summarized  in  the  following  proposition. 


Proposition  1  Consider  a  mechanical  system  with  equations  of  motion  prescribed  by 

-0 

dt  dq  dq 


where  C  is  given  by  (10).  Let  Q  =  F{q^x)  be  an  input  dependent  change  of  configuration  coordinates 
as  described  above.  Then,  in  Q -coordinates  the  equations  of  motion  take  the  same  form: 


dt  dQ  dQ 


(17) 


where  C  is  given  by  (13),  with  M,  A,  and  V  given  by  (I4),  (15),  and  (16)  respectively. 


Proof  It  is  straightforward  to  show  that  the  reduced  Lagrangian  has  the  form  claimed  when  ex¬ 
pressed  in  terms  of  Q-coordinates.  It  is  a  slightly  tedious  but  also  straightforward  calculation  to 
show  that  the  equations  of  motion  are  given  by  (17).  Alternatively,  if  the  reduced  Lagrangian  (11) 
arises  through  the  reduction  process  we  have  described,  it  may  be  shown  that  the  claimed  invariance 
follows  from  the  invariance  of  the  classical  Lagrange’s  equations.D 

3  Flat  Systems  and  Systems  with  Flat  Inputs 

This  section  proposes  some  basic  elements  of  a  theory  of  normal  forms  for  the  velocity  and  accel¬ 
eration  controlled  mechanical  systems  which  were  introduced  in  the  previous  section.  We  start  by 
considering  an  uncontrolled  mechanical  system  which  is  simple  in  the  sense  that  it  has  a  Lagrangian 
of  the  form 

L  =  2  9^Af(9)9  -  V (q),  (18) 

comprised  of  a  kinetic  and  potential  energy  term.  (Cf.  Abraham  and  Marsden,  1978.)  The  results  we 
present  are  steps  in  a  program  aimed  at  developing  an  intrinsically  second  order  theory  of  mechanical 
control  systems. 

In  the  study  of  first  order  nonlinear  control  systems  of  the  form  (1),  important  invariant  quali¬ 
tative  features  of  the  dynamics  are  revealed  by  computing  Lie  brackets  of  vectorfields  in  the  family 
{/(•jx)  :  X  €  U}.  We  refer  to  Nijmeijer  and  van  der  Schaft  [1990]  for  basic  results  on  Lie  bracket 
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conditions  related  to  controllability,  accessibility,  feedback  linearizability,  etc.  The  theory  of  first 
order  control  systems  also  applies  to  the  second  order  systems  treated  in  the  previous  section  if  we 
represent  them  in  the  usual  way  as  first  order  systems  (of  twice  the  dimension)  with  generalized 
forces  playing  the  role  of  control  inputs.  Our  goal  in  the  present  section  is  to  describe  an  intrinsi¬ 
cally  second  order  approach  in  which  curvature  provides  an  invariant  measure  of  nonintegrability. 
It  is  worth  remarking  that  while  there  has  been  a  great  deal  written  over  the  past  25  years  making 
connections  between  analytical  mechanics  and  nonlinear  control,  a  question  which  remains  almost 
completely  unanswered  is  whether  there  can  be  found  an  intrinsically  second  order  theory  of  con¬ 
trolled  mechanical  systems.  Almost  the  last  word  on  this  subject  in  fact  appeared  in  Brockett’s 
NASA  Ames  lectures. 

The  section  discusses  the  role  of  curvature  in  three  different  settings.  First,  we  recall  the  basic 
result  that  the  vanishing  of  the  curvature  tensor  associated  with  the  Levi-Civita  connection  of 
an  uncontrolled  simple  mechanical  system  (18)  is  necessary  and  sufficient  for  the  existence  of  a 
coordinate  system  with  respect  to  which  the  inertia  tensor  is  the  identity  matrix.  Second,  for  the 
case  of  a  simple  mechanical  system  in  which  every  degree  of  freedom  is  independently  actuated, 
we  recall  that  in  the  Hamiltonian  setting  the  same  curvature  vanishing  implies  the  existence  of  a 
(canonical)  coordinate  transformation  together  with  a  state  feedback  transformation,  not  involving 
conjugate  momentum  variables,  transforming  the  dynamics  into  a  decoupled  set  of  second  order 
integrators.  Finally,  we  consider  generalized  Lagrangian  control  systems  of  the  form  introduced  in 
the  previous  section.  For  these  systems,  there  is  defined  an  input  connection  and  related  curvature¬ 
like  quantities.  The  vanishing  of  both  the  curvature  defined  by  the  Levi-Civita  connection  and 
the  curvature  defined  by  the  input  connection  is  necessary  for  the  existence  of  an  input-dependent 
change  of  configuration  coordinates  under  which  the  inertia  tensor  is  transformed  to  the  identity 
matrix  and  all  terms  depending  on  input  accelerations  in  the  dynamics  are  eliminated.  Thus  for 
the  subclass  of  generalized  Lagrangian  control  systems  for  which  the  Levi-Civita  connection  has 
zero  curvature,  the  vanishing  of  the  curvature  of  the  input  connection  is  a  necessary  condition  for 
an  acceleration  controlled  mechanical  system  to  be  reducible  through  a  change  of  coordinates  to  a 
velocity  controlled  mechanical  system. 

The  simple  Lagrangian  (18)  gives  rise  to  the  Lagrangian  dynamics: 

n  n  51/ 

^  ^  ^  ^  +  ~Q~  =0,  (i  =  l,...,7l),  (19) 

where  the  Tkji  are  the  Christoffel  symbols  of  the  first  kind  introduced  in  the  previous  section.  The 
terms  they  define  in  equation  (19)  describe  the  inertial  forces  (Coriolis  and  centrifugal)  which  affect 
the  system.  It  is  clear  from  this  equation  that  if  M  does  not  depend  on  q,  Tkji  =  0  for  all  k,j,  i,  and 
these  inertial  forces  will  not  play  a  role  in  the  system’s  dynamics.  An  interesting  question  is  “When 
can  we  find  a  change  of  coordinates  such  that  in  the  new  coordinate  system  the  inertia  matrix  is 
constant?”  Note  that  the  inertia  matrix  M  is  constant  (i.e.  does  not  depend  on  the  configuration 
variables  q)  if  and  only  if  there  is  a  system  of  generalized  coordinates  (yi, . . .  ,yn)  in  terms  of  which 
the  kinetic  energy  is  expressed  as  |  SILi  Vi-  (Froo/.*  In  the  case  that  M  is  constant,  we  let  Y 
denote  the  constant,  symmetric,  positive  definite  matrix  square  root  of  M  and  let  y  =  Yq.  Then 
y  =  Yq,  and  the  kinetic  energy  is  reverse  implication  requires  no 

proof.) 

Hence,  whether  we  can  find  the  desired  change  of  coordinates  amounts  to  whether  there  exists 
a  diffeomorphism  y  =  F{q)  such  that  ^q'^Mq  =  ||ly|P*  Suppose  such  an  F  exists.  Then  y  = 

and  =  M.  For  such  an  M  we  compute 

dmij  ^  A  .  d^Ft  dFt  dFt  d^Fi  . 

^Qk  ^^^dqidqjdqj  dqi  dqjdqj^ 
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From  this  it  follows  from  an  easy  calculation  that 

^kji  = 


As  defined  in  Section  2,  the  Christoffel  symbols  (of  the  first  kind)  Tkji  are  associated  with  a  con¬ 
nection  (the  so-called  Levi-Civita  connection  defined  by  the  inertia  tensor  M)  and  a  Riemannian 
curvature  tensor.  It  is  a  classical  result  that  this  tensor  vanishes  if  and  only  if  there  is  a  diffeomorphic 
change  of  coordinates  Q  =  F{q)  such  that  the  inertia  tensor  may  be  factored  as  M  =  Jp  ■  Jp  where 
Jf  is  the  Jacobian  of  F.  To  see  this,  we  first  write  down  the  corresponding  Christoffel  symbols  of 
the  second  kind: 

i=l 

The  Riemannian  curvature  tensor  is  prescribed  by  the  corresponding  Riemann  symbols  of  the  first 
kind:  ^ 

=  ^  -  ^  +  E  (rj,r,.,  -  .  (21) 

If  F{-)  is  a  diffeomorphism  (local  coordinate  transformation)  such  that  M  —  Jp  ■  Jp, 

dVjti  ^  A.  d^F,  dF,  d^F,  d^F,  . 

9qk  ^ydqkdqjdqe  dqi  dqjdqt  dqidqj ' 


K 

n 

E 


9mij  druik 
dqk  dqj 
d^Fe  dFt 
^  dqjdqk  dqi ' 


dmkj 

dq, 


(20) 


and  a  permutation  of  indices  gives  a  similar  expression  for  .  It  is  straightforward,  although  a 
little  tedious,  to  show  that 


n 


n 


<7  =  1  )0=1 


d'^Fp  d'^Fp 
dqjdqk  dqidqt 


d^Fp  d^Fp  ' 
dqjdqe  dqidqk 


Hence  Rijki  =  0,  proving  that  the  curvature  vanishing  is  a  necessary  condition  for  there  to  be  a 
factorization  of  the  desired  type.  For  a  proof  of  sufficiency  and  a  broader  discussion  of  curvature 
and  its  vanishing,  the  reader  is  referred  to  Kobayashi  and  Nomizu  [1963]  or  Wolf  [1972]. 


Remark  3  The  Riemannian  curvature  may  equivalently  be  characterized  in  terms  of  Riemann 
symbols  of  the  second  kind: 


= 


It  is  not  difficult  to  prove  the  very 
second  kinds: 


direct  relationship  between  Riemann  symbols  of  the  first  and 
R/k(  =  Ern^^R<rike. 


Remark  4  While  there  are  components  of  the  curvature  tensor,  these  are  not  all  independent, 
and  the  Riemann  symbols  associated  with  the  Levi-Civita  connection  can  be  shown  to  satisfy: 


.  Rijke  +  Riktj  +  RiCjk  —  Oj 


10 


^ijkl  —  ^ijCk  j  S’lld 
Rijki  “  Rkiij' 

(See  Dubrovin  et  al  [1984],  Theorem  30.2.1.)  It  is  not  difficult  to  show  that  these  symmetries 
imply  that  only  n^(n^  -  1)/12  of  the  Riemann  symbols  of  the  first  kind  are  independent. 

Definition  4  A  simple  mechanical  system  of  the  form  (19)  will  be  called  fiat  if  the  corresponding 
Riemannian  curvature  tensor  (i.e.  the  set  of  Riemann  symbols  of  the  first  kind)  is  zero. 


Bedrosian  [1992]  and  Spong  [1998]  have  noted  that  this  type  of  flatness  points  to  simplified 
control  designs  in  the  case  of  fully  actuated  Lagrangian  and  Hamiltonian  systems.  We  review 
the  basic  ideas  and  then  turn  our  attention  to  the  more  complex  (and  more  interesting)  case  of 
super-articulated  (or  underactuated)  systems.  Consider  a  mechanical  control  system  specified  by  a 
Lagrangian  L{q,q)  =  ^q^M{q)q  -  V{q)  and  equations  of  motion 


dt  dqi  dqi 


m  {i  —  1, . . . ,  7i), 


where  each  Ui(*)  is  a  piecewise  continuous  control  input  (function  of  time  defined  on  a  suitable 
interval).  The  conjugate  momentum  vector  for  this  system  is  p  =  M{q)q  and  in  terms  of  the 
Hamiltonian  defined  by  the  Legendre  transformation 


H{q,p)  =  fp-L{q,q) 

=  ^P^M{q)-^p  +  V{q), 


(22) 


we  define  the  corresponding  Hamiltonian  control  system 


Q 

P 


dp 

dH 


(23) 


It  has  been  noted  (See,  e.g.  Bedrosian  [1992]  and  Spong  [1998].)  that  if  a  system  is  flat  there  is  a 
canonical  transformation  together  with  a  control  and  feedback  transformation  such  that  in  terms 
of  the  new  phase  space  and  control  variables,  the  system  (23)  takes  the  form  of  a  system  of  double 
integrators. 

To  understand  this,  suppose  the  mechanical  system  is  flat.  Let  M^{q)  be  the  n  x  n  positive 
definite  symmetric  square  root  of  the  inertia  tensor  M{q).  If  there  exists  a  function  F{q)  such  that 
^  =  M5(g)  as  above,  then  the  coordinate  transformation  Q  =  F{q)y  P  =  M^{q)q  is  canonical — in 
the  sense  that  it  preserves  the  symplectic  form.  We  refer  the  reader  to  Marsden  and  Ratiu  [1994]  or 
Arnold  [1989]  for  a  discussion  of  canonical  transformations.  That  this  transformation  is  canonical  is 
easily  seem  if  we  write 

P-dQ  =  q^M^  (q)  ■  ^dq 

=  q’^M{q)dq 
=  p-dq. 


We  rewrite  the  Hamiltonian  in  terms  of  Q,P“COordinates: 

Hiq,p)  =  ^p'^M{q)-^p  +  Viq) 

=  ^fMHq)MHq)q  +  V{q) 

=  l\\Pf  +  V{Q), 

■  =  n{Q,P), 
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where  V  is  defined  by  V{Q)  —  V{q). 


Proposition  2  In  terms  of  the  variables  Q,P,  the  equations  (23)  may  be  written  as 

Q  =  P 

P  =  +  (24) 

Under  the  feedback  u  =  M^v  +  Ms  which  is  independent  of  the  conjugate  momentum  variables 
P,  this  Hamiltonian  control  system  is  transformed  into  the  system  of  double  integrators 


Q  =  P 

P  =  V  (25) 


Proof:  The  coordinate  change  associated  with  the  Legendre  transformation  may  be  written 
explicitly  as 


Q  = 


P  = 


(26) 

(27) 


(since  We  wish  to  find  the  equations  of  motion  in  terms  of  Q,P-coordinates. 

Differentiating  (26),  it  is  clear  that  Q  =  P.  Differentiating  (27),  we  find  that 


P 


1-f— )' 

^dt  \dq  y 


']p+{ 


^N-l 
dq) 
dF 


P 

rd/dF\-K  rdF\-idH  rdF\ 

^dt  \da)  \  So  )  So  V  So  7 


rSF\-i 


^dt  \dq  J 


dq  J  dq  \dq  J 


u. 


Using  the  definition  of  the  coordinate  transformation  (26)- (27)  together  with  the  symmetry  proper¬ 
ties  of  ^  =  M2(g),  it  can  be  shown  that 


Prom  this  it  follows  that 

^dt\dq)  ^  \dq)  dq  dQ  dQ' 

proving  the  first  part  of  the  proposition.  The  second  part  is  clear.  □ 

More  generally,  a  system  of  the  form  (7)  or  (12)  is  said  to  be  fiat  if  the  Riemannian  curvature 
tensor  prescribed  by  the  Christoffel  symbols  vanishes.  In  terms  of  the  Lagrangian  system  (10), 
if  the  inertia  tensor  M  does  not  depend  on  the  input  x,  then  the  vanishing  of  the  curvature  tensor 
is  necessary  and  sufficient  for  there  to  be  a  change  of  generalized  coordinates  q  =  F{q)  such  that 
the  inertia  tensor  expressed  in  terms  of  the  g— coordinates  has  entries 


r  1 
1  0 


For  the  case  in  which  M  =  >W(g,2;),  explicitly  depends  on  the  input  x,  it  is  still  possible  to  have 
the  Riemannian  curvature  tensor  vanish.  (Of  course,  in  this  case,  the  Riemann  symbols  themselves 
will  have  no  dependence  on  the  variable  x.)  This  vanishing  is  a  necessary  and  sufficient  condition 
for  there  to  be  an  x-dependent  coordinate  transformation  with  respect  to  which  M  is  transformed 
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into  the  identity  matrix.  Let  q  =  F{q,  x)  be  such  a  coordinate  transformation,  and  let  q  =  G{q,  x) 
be  its  i-dependent  inverse.  Then  M(q,x)  =  and  in  g-coordinates,  the  Lagrangian  (10) 

takes  the  form 

C{q,’q\x,x)  =  -\\h\?  ■\rX^A{q,xyq-V{q-,x,x),  (28) 


where 


and 


BO  BF'^ 

A{q,x)  =  A{G{q,x),x)—  -  —  , 


.dG  i.T-dF^dF 


V(g;  X,  x)  =  V{G{q,  i);  x,  x)  -  x^  .4—  - -x 


X, 


with  all  component  entries  in  being  evaluated  at  g  =  G{q,x).  An  interesting  question  is  “When 
do  the  coupling  terms  A{q,x)  vanish?”  The  following  definition  is  useful  in  formulating  the  answer. 


Definition  ^  (i)  For  the  second-order  generalized  control  system  (12),  we  say  that  the  hatted 
symbols  define  the  input  connection.  The  input  connection  is  said  to  be  flat  if  the  associated 
Riemann  symbols  of  the  second  kind 


Rf 


li 


dxi 


Qp7  n 
fc=i 


dxi 


vanish  for  all  i,  £  =  1, . . . ,  m;  j,  7  =  1, . . . ,  n.  (ii)  Given  the  second  order  generalized  control  system 
(12),  we  say  that  the  pair  {x,x)  constitutes  a  set  of  flat  inputs  if  the  input  connection  is  flat.D 

The  desired  result  on  the  structure  of  the  Lagrangian  control  system  (12)  is  the  following. 


Theorem  1  Consider  the  Lagrangian  control  system  (10).  Let  U  xV  C  IR"  x  IR"*,  and  suppose 
that  F  :U  xV  -¥  IR"  is  an  input-dependent  change  of  coordinates  g  =  F{q,  x)  such  that  for  each  x 
the  metric  tensor  M  expressed  in  g-coordinates  has  Sij  as  its  ij-th  entry.  Suppose,  moreover,  that 
in  the  g-coordinates  all  cross  coupling  terms  .4(g,a;)  vanish.  Then  the  system  (defined  by  (10)) 
is  flat  and  has  flat  inputs,  which  is  to  say  both  the  hatted  and  unhatted  Riemann  symbols  of  the 
second  kind. 


Br'^  Br^  ^ 

+  Pk  -  A  Pi'-  ak) 

=  l,...,n. 


and 


R 


BV'^  BV^  ^ 


vanish. 

The  proof  of  this  theorem  uses  the  following  lemma. 

Lemma  1  Let  U  xV  C  IR”^  x  H”  be  some  neighborhood  of  (0,0),  and  let  fjiUxV-^  IR”  be  a 
smooth  mapping  for  j  =  1, . . . ,  m.  Given  q  ^V,  there  is  defined  a  neighborhood  FF  of  0  G  IR*”  and 
a  smooth  function  ^  :  W  IR”  such  that 

(i)  9{0)  =  q,  and 
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(ii)  =  /;•('’> ffW)  for  all  r  6  ty,  j  =  1, . . .  ,m 

if  and  only  if  there  is  a  neighborhood  of  (0,  q)  on  which 


dfj_  __  ^  /it  _  ^  fk\  _  Q 


dn  dvj  '  dqk‘ 

This  lemma  is  proved  in  Spivak,  [1970].  We  can  now  prove  Theorem  1. 


Proof  (of  Theorem  1).  Let  q  =  F(g,x)  be  as  in  the  hypothesis  of  the  theorem.  Then  using 
Proposition  1  at  the  end  of  Section  2,  the  g-coordinate  rendering  of  the  Lagrangian  (10)  is 

Ciq,  'q-,  X,  v)  =  M  {q,  x)q  +  v'^A{q,  x)q-  V{q;  x,  v) 

where  the  various  terms  correspond  to  those  given  by  formulas  (14)-(16). 

The  hypothesis  that  M  has  coordinate  entries  Sij  is  equivalent  to  writing  ^  =  I,  and 

since  ^  we  conclude  that  M(q,x)  =  The  argument  used  at  the  beginning  of  this 

section  shows  that  the  Reimannian  curvature  tensor  associated  with  this  M{q,x)  is  therefore  zero. 
The  second  part  of  the  hypothesis  is  the  statement  that 

dG^ 

—  7Vl(gr,a:) +  .4(g,a:)  =  0. 

Using  the  fact  that 

dG.„,  .  -  dF-\  - 

^(F(,,x),x)  = -^  (g.x) 

together  with  the  assumed  factorization  M  =  some  straightforward  algebra  shows  the  above 

equation  is  equivalent  to 

A{G{q,x),x)-^{G{q,x),x)  =  ^  {G{q,x),x), 


which  is  in  turn  equivalent  to  writing 


Mn 

in  terms  of  the  original  coordinates.  The  componentwise  rendering  of  this  equation  is 


(29) 


ati 


=  E 

jfe=i 


dFk  dPfc  £  =  1, . . .  ,m 
dxi  dqi  2  =  1, . . .  ,n. 


Similarly,  under  the  hypothesis  of  the  theorem, 

A  dPk  dFk 


mu 


E 


^  9qi  dqj 


i,j  = 


It  is  a  straightforward  calculation  to  show  that 


drukj  dag 
dxi  dqj 


daij  _  2  Fg  dFg 

dqk  “  dqjdxi  dqk  ' 
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We  have  called  this  quantity  the  input  symbol  of  the  first  kind  tejk  and  we  define  the  associated 
input  symbol  of  the  second  kind  by 

A:=l 

where  is  the  fik-th  entry  in  Using  these  definitions,  one  can  show  that 


dFi 

d(lj 

\ 

n 

( 

OFx 

dQk 

\ 

: 

-  E 

1 

dqj 

1 

k=l 

\ 

dha. 

dQk 

y 

dxi 


This  is  a  partial  differential  equation,  conditions  for  the  solution  of  which  are  covered  by  Lemma 
1.  Specifically,  for  each  Fi  we  are  interested  in  conditions  under  which  there  exists  a  function 
g  :  IR^  satisfying 

where 

f}ix,z)  = 

7=1 

According  to  Lemma  1,  a  necessary  condition  for  the  solution  of  (30)  is  that 


(30) 


dfe  dfi  ■^,dft  k  dfi  j.,  _ 


This  is  rendered  (componentwise) 


)=0. 


7=1  '  -  ky- 

Since  this  must  hold  identically  in  z  =  (zi, . . .  ,2:n),  we  have  the  desired  conclusion  that 


k=l 


Remark  5  In  light  of  the  identity  noted  in  Remark  3,  it  is  clear  that  the  vanishing  of  both  the 
Riemann  Symbols  of  the  first  kind  and  input  symbols  of  the  first  kind  is  similarly  a  necessary 
condition  following  from  the  hypothesis  of  Theorem  1. 


Definition  6  If  there  is  a  choice  of  coordinates  q  such  that  the  inertia  tensor  expressed  in  q- 
coordinates  is  the  identity  and  the  cross  coupling  terms  A{qj  x)  vanish,  we  refer  to  the  ^-coordinates 
as  flat  coordinates. 

Remark  6  The  vanishing  of  all  the  Riemann  symbols  (of  the  first  or  second  kind)  may  be  shown 
to  be  sufficient  for  the  existence  of  an  input-dependent  change  of  coordinates  such  that  the  inertia 
tensor  M{q,x)  is  transformed  into  the  identity  matrix.  The  argument  can  be  constructed  along 
classical  lines,  and  the  reader  is  referred  to  Spivak  [1970]  for  details. 

When  the  curvature  associated  with  the  Levi-Civita  connection  of  the  reduced  inertia  tensor  M 
(equation  (21)  does  not  vanish,  it  is  not  presently  known  whether  the  vanishing  of  the  input  symbols 
is  sufficient  to  eliminate  the  appearance  of  input  accelerations  in  (12). 


15 


Remark  7  Sections  4  and  5  will  treat  the  connection  between  the  geometry  of  second  order  control 
systems  and  recent  results  on  averaging.  Many  (indeed  most)  of  the  hard  results  to  date  which 
describe  close  connections  between  the  qualitative  dynamics  of  averaged  and  nonautonomous  versions 
of  Hamiltonian  and  Lagrangian  control  systems  with  oscillatory  inputs  have  been  established  for  the 
case  of  single  input  systems.  We  note  that  the  input  connection  is  trivially  flat  for  single  input 
systems. 

To  develop  physical  intuition  regarding  the  vanishing  conditions  in  Theorem  1,  we  consider  some 
prototypical  controlled  mechanical  systems,  beginning  with  rotating  heavy  chains. 

Example  1  (Rotating  heavy  chains)  The  global  dynamics  of  a  heavy  (hanging)  chain  undergoing 
forced  rotations  about  a  vertical  axis  will  be  crucially  dependent  on  how  the  joints  between  the  links 
of  the  chain  constrain  the  possible  relative  motions  of  the  links.  This  dependence  has  been  explored 
in  Baillieul  and  Levi  [1991].  Further  insight  into  the  global  dynamical  effects  of  relative  motion 
constraints  may  be  developed  in  terms  of  the  curvatures  defined  above.  Consider  the  two  rotating 
pendulum  systems  in  Figure  1.  Both  pendula  undergo  controlled  rotations  about  the  vertical  axis. 
In  Figure  1(a),  the  joint  by  which  the  pendulum  is  suspended  is  a  single  degree-of-freedom  revolute 
joint,  while  in  Figure  1(b),  the  link  is  suspended  by  a  (two  degree-of-freedom)  universal  joint.  The 
respective  Lagrangians  are 

2 

^ {a?  CQS^  <j)  +  siv?  +mgC cos (f> 

2  2t 


and 


=  (12  +  ^(H- sin^^))(A2-t-(^-t-Z2)^2 


+  (i —  cos<^  sm{2‘tp)^9  +  (a^  +  2P)  sixKpipO 

cos^  4)  cos^  Ip 

+{^  +  1^)  (sin^  (p  +  cos  <p  sin^  ‘tp)j  0^ 

+  mg£  cos  0  cos  ip. 

The  corresponding  dynamics  of  the  controlled  equations  are 


dt 


-  (a^  cos^  (p  + sin^  0)0j  = 


(^^ -I- +  ~ 

and  (in  abbreviated  form) 


(31) 


d  dL2 
dt  de 


=  u, 


d  dL2  dL2  _  Q 

dt  90  d<p 


d  dL2  dL2 
dt  dip  dip 


=  0. 


(32) 


In  both  cases  a  torque  u  is  applied  to  control  the  angular  velocity  6  of  rotation  about  the  vertical 
axis.  (See  Baillieul  [1987]  and  Marsden  and  Scheurle  [1993]  for  more  information  about  rotating 
heavy  chains.)  Applying  our  formal  reduction  procedure,  (10)  is  rendered  respectively 


A(0.0;^)  =  y(^"  +  y)«A"-Vi(0;0) 


“^'12 


where 

and 


Vi  (0;  6)  =  -(a^  cos^  0  +  sin^  <p)9^  -  mgi  cos  0, 


(33) 


d^2{<p,4>,^,1p\0)  =  -(0,V')A42 


it) 


-f-  0  A2 


V2{<p,ip;e) 


(34) 


16 


where 


/  ;2™2^  ,  0 
M2  =  M2{<l>,i))  =  M  2  .  u 


and 


f  — cos (f)  cos'll)  sin 'ip 

V2{4>,i>\B)  =  -|(o2cos2  0cos2i/»+ (sin^(/>  +  cos^sin^V’))^^ 
-mp^  cos<^  cosV’- 


The  significant  difference  between  the  two  rotating  pendulum  systems  is  the  absence  or  presence 
of  the  coupling  terms  A  in  (33)  and  (34)  respectively.  The  rotating  planar  pendulum  (31)  feels 
the  influence  of  rotation  only  through  velocity  terms  involving  9^.  The  dynamics  of  the  universal 
joint  pendulum  (32),  however,  depend  on  both  the  angular  velocity  6  and  angular  acceleration  6.  A 
natural  question  which  arises  is  whether  there  is  a  change  of  coordinates  which  makes  the  dynamics 
(32)  dependent  only  on  9  and  not  on  the  acceleration  9.  Theorem  1  provides  a  partial  answer,  but 
before  considering  this  in  detail,  we  also  examine  the  Riemann  symbols  of  the  first  kind  for  the 
inertia  tensor  of  the  unreduced  rotating  universal  joint  system: 


AjiM 


A2{<l>,i>)  \ 
M2{<I>,^)  ) 


where  M2  and  A2  are  as  above,  and 


A/2  Ip)  =  cos^  (p  cos^  Ip  + 


^  (sin^  (p  +  cos^  (p  sin^  'ip)  . 


It  is  a  straightforward  calculation  (which  is  tedious,  if  you  don’t  use  computer  algebra)  to  show 
that  none  of  the  six  independent  Riemann  symbols  of  the  first  kind  associated  with  the  3  x  3  inertia 
tensor  M2{(p,'ip)  is  zero.  We  can  compare  this  with  the  results  obtained  by  evaluating  corresponding 
Riemann  symbols  and  inputs  symbols  for  the  reduced  Lagrangian  system  defined  by  M2{(p,'ip)  and 
A2{{(p,'ip)’  For  the  2x2  tensor  M2{(p,'ip),  the  curvature  tensor  is  defined  by  the  single  Riemann 
symbol 

The  input  connection  is  nonzero  in  that 

2  2 

f'112  =  (^^  +  ^)  cos  (^  -  (^  -  i^)  cos  (^(cos^  V'  -  sin^  ip)- 

It  is  obvious  from  the  definition,  however,  that  the  input  symbols  for  any  single  input  system,  such 
as  this,  are  going  to  be  zero. 

The  answer  to  the  question  of  whether  there  is  a  change  of  coordinates  which  eliminates  the 
dependence  of  the  equations  of  motion  on  input  acceleration  terms  is  not  completely  resolved  by 
these  calculations.  The  nonvanishing  of  the  curvatures  associated  with  both  the  3x3  inertia  tensor 
of  the  unreduced  system  and  the  2x2  inertia  tensor  of  the  reduced  system  imply  that  there  is  no 
change  of  coordinates  such  that  the  dynamical  dependence  on  input  accelerations  is  eliminated  and 
the  dynamics  of  the  configuration  variables  (p  and  'ip  are  also  decoupled.  Whether  one  can  simply 
eliminate  the  input  accelerations  by  a  choice  of  coordinates  in  which  there  is  inertial  coupling  of  the 
^  and  'ip  dynamics  remains  for  the  moment  an  open  question. 


^)sin^2'0  \ 

+  f  [l  +  sin^V^]')/ 
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Example  2  Controlled  planar  body  with  pendulum  attachment  Here  we  consider  a  pair  of  rigid 
bodies  which  are  connected  by  a  simple  frictionless  single  degree-of-freedom  hinge  as  illustrated  in 
Figure  2(a).  One  of  the  bodies  (the  larger  one  depicted  in  the  figure)  is  assumed  to  have  actuators 
allowing  its  motion  in  the  plane  to  be  controlled  to  follow  any  prescribed  smooth  path.  No  actuation 
is  applied  directly  to  the  pendulum,  and  hence  it  moves  entirely  under  the  infiuence  of  gravity  and 
motion  of  the  controlled  body.  Take  as  generalized  coordinates  for  the  system  x  (the  horizontal 
displacement  of  the  large  body),  y  (the  vertical  displacement  of  the  large  body),  and  9  (the  angular 
displacement  of  the  pendulum  from  the  vertical,  downward  pointing  configuration).  The  3x3  inertia 
matrix  M  in  this  case  is  given  by 

(nia nrib  0  mt^cosO 

0  ma  +  mh  mb£  sin  0 
rubicosO  TUbi  sin  6  h 

The  sub“blocks  of  the  inertia  matrix  M  which  are  of  interest  are  A  —  A{6)  =  {mb^  cos  6,mbi  sin  6) 
and  M  =  hy  where  rua  is  the  mass  of  the  controlled  body,  is  the  mass  of  the  pendulum,  h  is 
the  pendulum  inertia  about  the  hinge  point,  and  ^  is  the  distance  from  the  hinge  point  to  the  center 
of  mass  of  the  pendulum.  Thus  in  the  reduced  Lagrangian  control  system,  x  and  y  are  the  inputs, 
while  6  is  the  configuration  variable  whose  motion  is  governed  by  the  reduced  Lagrangian  dynamics. 

The  Christoffel  symbols  of  the  first  kind  associated  with  the  Levi-Civita  connection  of  the  unre¬ 
duced  inertia  tensor  M  are  all  zero  except  for  Fasi  =  —rubisinO  and  r332  =  mbicosO.  Nevertheless, 
the  Riemann  curvature  tensor  (21)  vanishes,  and  from  this  we  conclude  that  there  is  a  choice  of 
coordinates  with  respect  to  which  the  inertia  tensor  takes  the  form  of  the  identity  matrix.  It  is 
also  easy  to  verify  that  the  vanishing  conditions  of  Theorem  1  are  satisfied.  Indeed,  the  reduced 
inertia  tensor,  Ad  =  /&,  is  a  scalar,  and  hence  the  Riemannian  curvature  is  trivially  zero.  The  input 
connection  tuj,  i,j  =  1,2  is  zero,  and  hence  the  input  curvature  is  also  zero. 

Example  3  Controlled  planar  body  coupled  with  extending  pendulum  attachment  The  concluding 
example  treats  another  system  for  which  the  curvature  vanishing  conditions  of  Theorem  1  are  not 
satisfied.  The  mechanism  depicted  in  Figure  2(b)  is  similar  to  the  previous  example,  with  the 
significant  difference  being  that  we  assume  there  is  an  internal  mechanism  which  causes  the  length 
of  the  pendulum  attachment  to  depend  on  the  x-coordinate  of  the  body.  Specifically,  we  assume 
the  length  of  the  pendulum  is  ax,  where  a  >  0  is  some  fixed  constant.  We  assume  again  that 
the  pendulum  is  attached  by  a  frictionless  hinge.  To  simplify  the  discussion  (and  with  no  loss  of 
generality)  we  idealize  the  model  so  that  the  pendulum  is  comprised  of  a  point  mass  rUb  and  a 
massless  linkage  between  the  point  mass  and  hinge.  The  inertia  matrix  M  in  this  case  is  given  by 

(  ra a mb TUbOt^ '\’2mbOt sin 9  — m^acos^  mbOLXCOs9 

-mba  cos  9  ma  +  mbax  sin  9 

mbOLX  cos  9  mbOLXsin9  mbOL^x*^ 

While  this  example  is  admittedly  somewhat  artificial,  it  serves  to  illustrate  some  of  the  features 
of  the  flatness  properties  we  have  discussed.  As  in  the  previous  example,  the  inputs  here  are  taken 
to  be  the  x  and  y  positions  of  the  controlled  body.  We  analyze  “nonintegrability”  in  the  dynamical 
relationship  between  the  input  variables  and  the  configuration  9  of  the  attachment  in  three  different 
ways.  First,  in  computing  the  components  of  the  Riemann  curvature  tensor  (21),  we  find  that  all 
six  independent  components  are  zero  except 

n  ^ _ m. _ 

1313  8  ((1  ^  2  a)  cos^  9  sin  9  mb'^  -H  (1  +  H-  sin  9)  mb  ma  H-  ma^) 

where 

N{9)  =  --((-1  +  2 a)  mb  (-2  (-l  +  2a)  cos^^m^^- 
(-3  +  cos(2  0)  -  4  sin0  “  4a^  sin9)  m^ma 
+  4  sin 
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This  means  of  course  there  is  no  change  of  (x,  j/,  0)-coordinates  such  that  the  unreduced  inertia 
tensor  M  has  the  form  of  the  identity  matrix.  Secondly,  we  analyze  the  reduced  Lagrangian  in  light 
of  Theorem  1.  Since  in  this  case  the  reduced  inertia  tensor  M{9;x,y)  =  is  of  dimension 

1  X  1,  we  shall  only  need  to  analyze  the  input  connection  to  verify  the  hypotheses  of  the  theorem. 
The  input  connection  has  two  components:  =  |  and  =  0.  That  is  unbounded  as  a:  0 

reflects  the  fact  that  the  inertia  tensor  is  singular  in  the  a;  =  0  limit.  Although  the  input  connection 
is  nonzero,  the  vanishing  of  all  the  relevant  curvatures  shows  that  it  is  flat  in  the  dense  of  Definition 
4.  Thus  the  necessary  conditions  of  the  theorem  are  satisfied.  Since  we  know  there  is  no  choice 
of  coordinates  in  terms  of  which  M  is  rendered  as  the  identity  matrix,  the  necessary  conditions  of 
Theorem  1  thus  fail  to  be  sufficient.  The  theory  at  present  does  not  answer  the  question  of  whether 
or  not  the  system  can  be  rewritten  in  such  a  way  that  the  dynamics  do  not  depend  on  accelerations 
of  the  inputs. 

There  is  a  third  and  final  sense  in  which  the  influence  of  the  inputs  {x,y)  in  this  example  is 
nonintegrable,  and  this  will  be  discussed  in  the  following  section. 

4  Averaging  Lagrangian  and  Hamiltonian  Systems  with  Os¬ 
cillatory  Inputs 

For  systems  such  as  Example  1(a)  where  we  may  express  the  Lagrangian  as  the  sum  of  a  simple 
kinetic  energy  term  plus  a  time-dependent  potential,  it  is  natural  to  conjecture  that  given  small- 
amplitude  high-frequency  periodic  inputs  x(-),  their  influence  is  approximately  equal  to  the  effect 
of  a  potential  force  determined  by  the  time-averaged  potential.  There  is  a  growing  body  of  theory 
concerning  the  conditions  under  which  this  conjecture  is  true.  (See  Weibel  [1997]  and  Weibel  and 
Baillieul  [1998]  for  details.)  The  role  of  flatness  and  flat  inputs  in  determining  the  dynamic  response 
of  (11)  undergoing  high  frequency  oscillatory  forcing  is  not  completely  understood,  but  it  is  clear 
that  the  analysis  can  be  simplified  considerably  for  systems  which  are  flat  in  the  sense  described. 
One  specific  simplification  which  occurs  for  the  case  of  flat  systems  with  flat  inputs  is  that  the 
averaged  potential  may  be  computed  (with  respect  to  flat  coordinates)  by  taking  simple  averages 
of  the  amended  potential  V{q]Xjv)  in  (10).  For  systems  for  which  we  cannot  find  flat  coordinates, 
there  seems  to  be  additional  complexity  and  additional  steps  in  the  analysis.  One  approach  which 
works  (at  least  formally)  in  the  general  case  is  to  move  to  the  Hamiltonian  setting.  Here  it  is  always 
possible  to  write  down  the  averaged  potential 

As  a  general  starting  point,  we  return  to  the  Lagrangian  control  system  (10)  which  we  rewrite 
here:  ^ 

C{q,q;x,v)  =  -q^M{q,x)q  +  v'^A{q,x)q  -V{q-,x,v). 

To  transform  to  a  Hamiltonian  formulation,  we  perform  a  Legendre  transformation,  writing  T-L  = 
p^q  -  C  and  expressing  Ti  in  terms  of  (g,p).  This  Ti  is  a  noncanonical  Hamiltonian  which  takes  the 
form 


'H{qyP\x,v)  =  ^{p- A^v)'^M  ^{p-A'^v)  +  V  (35) 

We  expand  equation  (35)  and  apply  simple  averaging  to  yield  the  averaged  Hamiltonian: 

?{(g,p)  =  ^p'^ M~^p  “  v'^AM~^p  +  ^  AM~^A^v  +  V.  (36) 

z  z 

Here  the  overbars  indicate  simple  averages  over  one  period  (T)  have  been  taken:  given  any  piecewise 
continuous  function  F{q,pfl),  the  simple  average  is  given  by  F  =  /q  F{q,pfl)dt,  where  for  the 
purpose  of  evaluating  this  integral  q  and  p  are  regarded  as  constants. 

For  H,  there  is  an  obvious  decomposition  into  kinetic  and  potential  energy  terms  in  the  case  that 
v'^AM-^  =  0: 
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n{q,p) 


(37) 


=  ^p^M~^p+^v'^AM~^A'^v  +  V. 

Zt  z 

' - V - ^  ^ - - - ' 

avg.  kin.  averaged  potential 

energy 

There  are  a  number  of  conditions  under  which  AM~^  =  0.  This  will  trivially  be  the  case  when 
we  are  able  to  find  a  choice  of  coordinates  with  respect  to  which  the  A  terms  do  not  appear  in 
(10).  In  Weibel  and  Baillieul  [1998]  it  was  shown  that  for  zero-mean  input  v(-)  exactness  of  the 
1-form  M~^(gix)A(q}X)  dx  is  also  sufficient  to  guarantee  that  v'^AM^^  =  0.  We  remark  that 
this  exactness  condition  is  satisfied  for  the  systems  in  Examples  1  and  2  but  not  Example  3  in  the 
previous  section.  _ 

In  the  case  that  v’^AM^^  7^  0,  it  remains  possible  to  decompose  the  averaged  Hamiltonian  (36) 
into  the  sum  of  averaged  kinetic  and  potential  energies,  although  the  description  becomes  more 
involved.  If  u  7^  0,  then  the  corresponding  input  variable  x{t)  will  not  be  periodic.  There  will 
be  a  “drift”  in  the  value  of  x{t)  which  changes  by  an  amount  v  •  T  every  T-units  of  time.  We 
may  nevertheless  study  averaged  Hamiltonian  systems  in  this  context.  We  rewrite  the  averaged 
Hamiltonian  (36)  as 

T-L{q^p)  =  AM~^p AM'~^  M~^A^v 

+  ^v'^AM-^Av-  ^v'^AM-^  M-^A'^v  +  V 

=  I  (j^P  -  (M^P  -  M-^A^v)  (38) 

> - - / 

averaged  kinetic  energy 

+  ^yT'AM-^AT'v  -  ^vT'AM-'^  M-^A^v  +  V. 

N - ' 

averaged  potential 

The  formal  distinction  which  appears  between  the  zero-mean  and  non-zero-mean  cases  ((38)  and 
(39)  respectively)  is  that  in  the  latter  case  both  the  averaged  kinetic  and  potential  energy  terms  are 
adjusted  to  refiect  the  net  (average)  motions  of  the  input  variables  (a;(-),u(*)).  There  are  important 
relationships  which  can  be  established  between  the  dynamics  associated  with  the  averaged  Hamilto¬ 
nian  (39)  and  the  dynamics  of  the  periodically  forced  system.  The  reader  is  referred  to  Weibel,  1997, 
Weibel  et  al,  1997,  and  Weibel  k  Baillieul,  1998  for  details.  It  is  important  to  mention  that  for 
mechanical  systems  in  which  u  7^  0  and  M  depends  explicitly  on  x  in  (35),  the  averaging  analysis  of 
this  chapter  may  not  provide  an  adequate  description  of  the  dynamics.  Indeed,  in  this  case,  ||a:(t)|| 
will  not  remain  bounded  as  t  00,  and  if  M{x{t),q2)  also  fails  to  remain  bounded,  the  averaged 
potential  will  inherit  a  dependence  on  time  which  will  make  it  difficult  to  apply  the  critical  point 
analysis  proposed  below.  Despite  this  cautionary  remark,  we  shall  indicate  how  our  methods  may 
be  applied  in  many  instances  where  v  ^0. 


5  stability  and  Flatness  in  Mechanical  Systems  with  Oscil¬ 
latory  Inputs 

The  Lagrangian  (10)  gives  rise  to  the  Lagrangian  dynamics  (11),  which  we  rewrite  here  as 

n  m  n  n  m 

^  ^  '^kjqj  ^  ^  ^ijkqiqj  "i”  ^  ^  ^  ^  ^Cjk'^iqj  “  ■^(^)j  (^  “  I5  •  •  •  >  ^)j  (^^) 

j=l  £=1  ij~l  i=l  i=l 
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where  Fijk  and  ttjk  are  defined  as  in  Section  2  in  terms  of  the  entries  and  rriij  in  the  m  x  n  and 
nxn  matrices  A{q,x)  and  M{q,x)  respectively.  We  seek  to  understand  the  stability  of  responses 
to  oscillatory  forcing  of  (11)  in  terms  of  the  corresponding  averaged  potential  which  was  introduced 
in  the  previous  section. 

VA{q)  =  M-^A^v  +  V.  (40) 

An  answer  to  this  question  is  provided  by  the  averaging  principle  for  Lagrangian  systems: 
Averaging  Principle  for  Periodically  Forced  Lagrangian  Systems:  The  dynamics  of  the 
periodically  forced  system  (11)  are  locally  determined  in  neighborhoods  of  critical  points  of  the 
averaged  potential  Va{q)  a-s  follows: 

1.  If  q*  is  a  strict  local  minimum  of  Vyi(-),  then  provided  the  frequency  of  the  periodic  forcing 
?;(•)  is  sufficiently  high,  (11)  will  execute  motions  confined  to  a  neighborhood  of  q*. 

2.  If  {q,p)  =  (^*,0)  is  a  hyperbolic  fixed  point  of  the  averaged  system  (11),  then  there  is  a 
corresponding  periodic  orbit  of  (11)  such  that  the  asymptotic  stability  properties  of  the  fixed 
point  (g*,0)  of  the  properly  averaged  version  of  (11)  coincide  with  the  asymptotic  stability 
properties  of  the  periodic  orbit  with  respect  to  (11). 

□ 

Roughly  but  simply  stated,  the  averaging  principle  for  forced  Lagrangian  systems  says  that  for 
periodic  forcing  functions  u(*)  of  sufficiently  high  frequency,  the  dynamics  of  the  averaged  system 
provide  a  good  local  approximation  of  the  dynamics  of  the  nonautonomous  (forced)  system.  Unfor¬ 
tunately,  this  averaging  principle  is  less  than  a  theorem,  and  it  has  been  rigorously  established  only 
in  special  cases.  We  highlight  several  of  these  in  the  present  section  and  refer  to  Baillieul  [1995]  and 
Weibel  and  Baillieul  [1998]  for  further  details.  The  main  idea,  of  course,  is  that  for  sufficiently  high 
frequencies  the  effect  of  forcing  (11)  with  an  oscillatory  input  u(-)  will  be  to  produce  stable  motions 
confined  to  neighborhoods  of  relative  minima  of  Va(*)-  A  naive  approach  to  stability  analysis  is  to 
linearize  the  dynamics  (11)  about  such  relative  minima.  Unfortunately,  it  is  only  in  special  cases 
that  such  linearizations  captures  the  stabilizing  effects  implied  by  an  analysis  of  the  averaged  poten¬ 
tial.  Indeed,  it  turns  out  that  in  general,  the  averaged  potential  depends  on  second  order  jets  of  the 
coefficient  functions  A{q)  and  M{q).  (We  continue  to  assume  that  M  and  A  depend  on  the  input 
X  in  general,  but  we  suppress  this  variable  to  simplify  notation  in  the  present  section.)  To  further 
simplify  the  presentation,  we  shall  restrict  our  attention  to  the  case  of  zero-mean  oscillatory  forcing 
in  which  M~^A^v  =  0. 

Suppose  go  is  a  strict  local  minimum  of  (40).  Applying  a  high-frequency  oscillatory  input  u(-),  we 
shall  look  for  stable  motions  of  (11)  in  neighborhoods  of  (g,  g)  =  (go,0).  Of  course,  even  when  there 
are  such  stable  motions,  (go,0)  need  not  be  a  rest  point  of  (11)  for  any  choice  of  forcing  function 
u(*).  (Cf.  the  “cart-pendulum”  dynamics  treated  in  Weibel  et  ah  [1997].) 

To  analyze  the  relationship  between  (11)  and  (40)  in  more  detail,  let  us  assume  go  =  0.  This 
assumption  is  made  without  loss  of  generality,  since  we  may  always  change  coordinates  to  make  it 
true.  Write 

A{q)  =  Ao  4*  Ai{q)  +  A2{q)  +  h.o.t.; 

M{q)  =  Mo  +  Mi{q)  H”  M2(q)  h.o.t., 

where  the  entries  in  the  nxn  matrix  Mkio)  are  homogeneous  polynomials  of  degree  k  in  the 
components  of  the  vector  g,  and  similarly  for  the  mxn  matrix  ^^(g).  The  averaged  potential  has 
a  similar  expansion,  which  we  write: 

V^(g)  =  Vo  +  Vi(g)  +  V2(g)  +  h.o.t., 

All  terms  in  all  these  expressions  may  depend  on  the  input  x,  but  we  do  not  explicitly  display  this 
dependence.  We  note,  however,  that  each  term  in  the  expansion  of  Va  is  the  sum  of  a  term  of  the 
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same  order  in  the  original  potential  plus  an  algebraic  function  of  entries  in  terms  up  to  the  given 
order  in  the  expansions  of  A  and  M. 

There  are  two  essentially  different  cases  which  can  be  described  in  terms  of  these  expansions. 
The  case  that  ^  0  is  qualitatively  quite  different  from  the  case  in  which  =  0.  It  is  associated 
with  qualitative  dynamical  behavior  we  have  referred  to  as  hovering  motions.  These  are  motions 
confined  to  neighborhoods  of  critical  points  of  the  averaged  potential  which  are  not  equilibrium 
points  of  the  corresponding  forced  or  nonautonomous  system.  We  will  not  treat  this  case  further, 
and  refer  to  Weibel  and  Baillieul  [1998]  for  further  details.  We  suppose  in  the  remainder  of  the 
section  that  Aq  =  0,  in  which  case  0  is  both  a  critical  point  of  the  averaged  potential  and  a  rest 
point  of  the  forced  system  (11).  We’ll  show  that  in  this  case  the  averaged  potential  depends  only  on 
first  order  jets  of  the  coefficients  of  (11)  when  A{qo)  =  0.  In  Weibel  and  Baillieul  [1998]  it  has  also 
been  shown  that  the  condition  -4(^0)  =  0  is  also  necessary  and  sufficient  for  the  local  minimum  qo 
of  the  averaged  potential  to  define  a  corresponding  fixed  point  (rather  than  a  periodic  orbit)  of  the 
forced  (nonautonomous)  Hamiltonian  system  associated  with  (10). 

Slightly  refining  our  notation,  let  A^{q)  denote  the  ^-th  column  of  the  n  x  m  matrix  A'^{q).  Then 
we  have 

A^{q)  =  A{  '  q  +  (terms  of  order  >  2),  and 
Miq)  —  Mo  +  (terms  of  order  >  1), 

where  we  interpret  Mq,  ,  A'^  as  n  x  n  coefficient  matrices.  Using  the  same  notation,  we  also 

have  an  expansion  of  the  (original)  potential  function: 

Viq)  =  ^q^V2q  +  h.o.t. 


We  recall  the  following  result. 

Proposition  3  Suppose  u(*)  is  an  -valued  piecewise  continuous  periodic  function  of  period  T  >  0 
such  that  V  =  ^  fo  v{s)ds  —  0.  Suppose,  moreover,  that  Aq  —  0.  Then  the  averaged  potential  of 
the  Lagrangian  system  (11)  agrees  up  to  terms  of  order  2  with  the  averaged  potential  associated  with 
the  linear  Lagrangian  system 


m 

Moq  +  Yl  +  V2  •  9  =  0.  (41) 

i=i 

Proof;  The  proof  follows  from  examining  a  series  expansion  of  Va  about  the  point  g  =  0.  We  refer 
the  reader  to  Weibel  and  Baillieul  [1998]  for  details. □ 

There  are  several  points  worth  noting  here.  First,  (41)  is  exactly  the  usual  (tangential)  lineariza¬ 
tion  of  (11)  with  respect  to  the  state  variables  (g,g)  about  the  rest  state  (g,  g)  =  (0,0).  Also,  if 
contrary  to  the  hypothesis  Ao  ^  0,  then  such  a  linearization  does  not  give  a  valid  local  approxima¬ 
tion  to  the  dynamics  of  (11),  and  the  Hessian  of  the  averaged  potential  (40)  involves  terms  of  higher 
order  in  both  A  and  M. 

A  deeper  connection  with  stability  is  now  expressed  in  terms  of  the  following  theorem. 

Theorem  2  Suppose  w{')  is  an  -valued  piecewise  continuous  periodic  function  of  period  T  >  0 
such  that  ^  /o^  w{s)ds  =  0.  Consider  the  linear  Lagrangian  system  (4I)  with  input  v{t)  = 
w{ujt),  and  suppose  A^^  =  A{  for  £  =  1, . . .  ,m.  The  averaged  potential  for  this  system  is  given  by 

-  m 

VA{q)  =  lq^[V2  +  Y.  (42) 

ij=l 

where  aij  =  (l/T)  Wi{s)wj{s)  ds.  If  the  matrix  is  positive  definite,  the  origin  (g,  g)  =  (0, 0) 
of  the  phase  space  is  stable  in  the  sense  of  Lyapunov  provided  u)  is  sufficiently  large. 
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□ 


This  theorem  has  been  proved  in  Baillieul  [1995].  Its  proof  apparently  requires  the  assumption 
that  the  coefficient  matrices  A{  are  symmetric.  Suppose  M(q,x)  =  M(q)  does  not  depend  on  x. 
Then 


and  f«j*;(0)  is  the  jA:-th  element  in  the  skew  symmetric  matrix  -  A{.  Thus,  in  this  case,  the 
hypothesis  that  A{  is  symmetric  for  £  =  1, ...  ,m  is  equivalent  to  assuming  the  input  connection 
vanishes  at  the  rest  point.  If  A{  is  not  symmetric,  the  vanishing  of  the  input  symbols  Rju  is 
equivalent  to  the  set  of  matrices 


{{Af  -  A[)Mq^  :i  =  l,...,Tn} 


being  commutative.  These  observations  are  interesting,  because  it  is  precisely  in  cases  where  flatness 
is  present  that  we  are  able  to  prove  the  stability  of  motions  based  on  an  analysis  of  the  averaged 
potential. 


6  Concluding  Remarks 

Like  many  of  the  chapters  in  this  volume,  the  present  one  treats  the  interface  between  nonlinear 
control  theory  and  mechanics.  Second  order  control  systems  are  obviously  central  in  this  regard,  and 
the  intrinsic  geometric  approach  we  have  presented  builds  naturally  on  earlier  work  of  Roger  Brock- 
ett.  We  have  shown  that  in  using  the  language  of  affine  connections,  it  is  possible  to  state  conditions 
involving  the  vanishing  of  curvature-like  quantities  that  must  be  satisfled  in  order  for  acceleration 
controlled  Lagrangian  systems  to  be  rewritten  (more  simply)  as  velocity  controlled  Lagrangian  sys¬ 
tems.  These  geometric  conditions  also  bear  on  the  simplicity  of  computing  the  averaged  potential 
and  assessing  the  stability  characteristics  of  such  systems  under  the  influence  of  high-frequency  os¬ 
cillatory  inputs.  It  is  of  interest  to  note  that  the  vanishing  of  the  input  connection  occurs  precisely 
in  the  case  that  we  have  been  able  to  prove  Lyapunov  stability  of  forced  “linear  Lagrangian”  sys¬ 
tems.  Since  our  vanishing  conditions  are  necessary  for  the  elimination  of  certain  Coriolis  terms  in 
the  unreduced  Lagrangian,  they  may  also  bear  on  questions  of  dissipation  induced  instability  which 
have  been  discussed  by  Bloch  et  al  [1994]. 
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Motor 


Figure  1:  (a)  Rotating  planar  pendulum  joint;  (b)  Rotating  universal  pendulum  joint.  We  consider 
a  simple  rotating  pendulum  suspended  by  either  type  of  joint  as  depicted.  The  pendulum  consists 
of  a  link  of  length  which  we  assume  to  be  massless  and  at  the  end  of  which  are  four  equal  masses 
distributed  as  shown  so  that  the  moment  of  inertia  of  the  pendulum  is  nonsingular. 
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Figure  2:  (a)  Controlled  planar  body  with  a  pendulum  attachment;  (b)  Controlled  planar  body  with 
an  extending  pendulum  attachment.  The  length  of  the  extensible  pendulum  is  ax^  where  x  is  the 
^-coordinate  of  the  hinge  point  connecting  the  two  bodies. 
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Abstract 

This  report  presents  new  results  on  the  behavior  of  an  n-link  pendulum  which  is  controlled  by 
means  of  an  oscillatory  input  forcing  one  end  link.  Drawing  on  recent  research  on  the  oscillatory 
control  of  second  order  nonlinear  systems,  we  adopt  the  averaged  potential  as  the  primary  tool 
in  our  equilibrium  and  stability  analysis.  While  such  control  of  single-degree^of-freedom  systems 
has  been  studied  in  fair  detail[15],  researchers  are  only  beginning  to  explore  the  issues  which 
arise  in  the  control  of  multi- degree-of-freedom  systems.  Of  specific  interest  is  the  behavior  when 
the  system  possesses  bifurcations  and  more  than  one  stable  equilibrium.  In  this  report,  we 
present  a  study  of  the  bifurcations  and  stability  of  a  periodically-forced  cart  and  n-pendulum 
on  an  inclined  plane.  Periodic  open-loop  forcing  of  the  single  degree-of-freedom  pendulum  has 
previously  been  shown  to  be  a  robust  and  effective  means  of  stabilizing  the  inverted  equilibrium, 
and  this  robustness  carries  over  to  the  multi-degree-of-freedom  case  provided  attention  is  paid  to 
matching  the  forcing  amplitude  to  the  characteristic  length  scales  of  the  system.  In  addition,  we 
present  numerical  results  on  the  basins  of  attraction  of  the  vertically  forced  double  pendulum. 
We  present  an  exact  model  for  the  n-degree-of-freedom  system,  nondimensionalize  the  model, 
and  compute  the  averaged  potential.  Equilibria  and  their  stability  chairacteristics  are  found 
through  a  bifurcation  analysis  of  the  averaged  potential.  For  a  double  pendulum  system,  we  are 
able  to  provide  a  detailed  comparison  of  analytical,  numerical,  and  experimental  results — all  of 
which  are  found  to  be  in  close  agreement. 


1  Introduction 

Recent  research  (e.g.  [10]  and  [6])  has  pointed  to  the  geometric  nature  of  stabilization  effects  pro¬ 
duced  by  oscillatory  forcing  of  controlled  mechanical  systems.  These  results  togehter  with  the 
nondimensionalized  theory  presented  in  [18]  suggest  that  high-frequency  oscillatory  control  designs 
can  be  effectively  applied  to  mechanical  systems  over  a  very  wide  range  of  characteristic  length 
scales.  In  particular,  there  is  great  appeal  in  applying  what  are  essentially  robust  open-loop  meth¬ 
ods  for  the  control  of  micromechanisms  where  they  allow  us  to  avoid  the  very  challenging  problem  of 
closing  feedback  loops  using  noisy  and  unreliable  sensor  technologies.  While  a  great  deal  of  previous 
work  on  the  oscillatory  control  of  mechanical  systems  has  been  aimed  at  so-called  superarticulated 
mechanical  systems  comprised  of  multibody  chains  (e.g.  [18], [17],  and  [15]),  recent  work  suggests 
interesting  applications  of  the  basic  concepts  to  the  control  of  micromechanical  systems  such  as 
bubble  dyanamics  ([9])  or  the  silicon  cantilevered  beam  and  plate  structures  arising  in  MEMS  ap¬ 
plications  ([7]).  A  crucial  result  in  applying  these  methods  to  systems  at  small  characteristic  length 
scales  is  that  stabilization  effects  are  scale  invariant.  The  present  report  is  aimed  at  illustrating 
some  of  the  fundamental  questions  which  are  raised  in  studying  this  issue,  and  we  have  chosen  to 
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to  the  Center  for  Dynamics  and  Control  of  Smart  Structures  (through  Harvard  University)  is  gratefully  acknowl¬ 
edged.  Partial  support  from  the  U.S.  Office  of  Scientific  Research  through  grant  F49620-96-1-0059  is  also  gratefully 
acknowledged. 
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place  the  discussion  in  the  context  of  the  response  of  heavy  chains  to  rapid  oscillatory  forcing.  This 
has  the  threefold  advantage  of  pushing  the  theory  further  into  the  real  of  multi-degree-of-freedom 
mechanical  systems,  giving  the  questions  of  length  scale  a  very  clear  physical  interpretation,  and  at 
the  same  time  making  contact  with  a  body  of  previously  published  work. 

Our  starting  point  is  to  consider  controlled  mechanical  systems  which  can  be  described  by  the 
Lagrangian 

C{Q,  Q;  v)  =  ^Q^M{Q)Q  +  vA{QfQ  -  Va(Q;  v),  (1) 

where  u  E  M  is  the  control  input,  M{Q)  E  is  a  state-dependent  inertia  tensor,  A{Q)  E  E” 
is  a  Coriolis  coupling  vector,  and  Va{Q\v)  E  E  is  known  as  the  augmented  potential  Such  systems 
typically  arise  from  n  +  1  degree-of-freedom  systems  where  the  velocity  of  one  degree  of  freedom 
is  viewed  as  the  control  input  v.  Details  of  the  model  order  reduction  process  leading  to  (1)  may 
be  found  in  [4].  The  noncanonical  Hamiltonian  'H{Q,P)  corresponding  to  (1)  is  found  through  the 
Legendre  transformation  P)  =  P^Q  -  C: 

H{Q,  P;  ^)  =  i(P  -  vA{Q)fMiQ)-\P  -  vA{Q))  +  Va(Q;  v)  (2) 

where  P  =  {Pi)P2,  •  •  • ,  Pn}^  ai'e  canonical  momenta  defined  by  P  =  dC/dQ.  We  refer  to  (2)  as 
“noncanonical”  because  for  an  arbitrary  input  u,  energy  is  not  necessarily  conserved  (i.e.  dH/dt  ^  0). 
If  we  assume  now  that  u  is  a  periodic  function  of  time  (i.e.  v{t  -h  T)  =  u(^)  for  some  T  >  0),  jn 
principle  HiQ^P’^v)  can  be  averaged  over  one  period  of  v{t)  to  obtain  the  averaged  Hamiltonian  Hy 
which  is  written 


niQ,  P)  =  5  (P  -  vA{Q)f  MiQ)-^  (P  -  vAiQ))  +  Va{Q),  (3) 

where  v  =  ^  v{t)dty  ^  v^{t)dty  and 

VAiQ)  A{Q)'^M{Q)-^A{Q)  +  VaiQ)  (4) 

is  known  as  the  averaged  potential  The  resulting  averaged  Hamiltonian  is  a  proper  Hamiltonian  (i.e. 
dH/dt  =  0),  and  the  equilibria  and  stability  of  the  associated  dynamics  may  be  studied  in  terms 
of  the  critical  points  of  Va{Q)-  The  similarity  of  the  dynamics  associated  with  (3)  to  the  dynamics 
associated  with  (2)  is  largely  dependent  on  the  forcing  amplitude  and  frequency.  At  the  moment, 
this  method  of  averaging  is  a  formal  method  which  has  been  shown  to  give  meaningful  results  for 
many  systems,  and  the  development  of  a  more  rigorous  theory  of  averaging  for  such  systems  is  an 
active  research  area. 

The  widely  studied  vertically-forced  pendulum[8,  4,  18]  possesses  a  Hamiltonian  of  the  form  (2), 
and  the  stabilization  of  a  more  general  pendulum  system  using  the  method  of  Hamiltonian  averaging 
has  recently  been  presented  in  [18].  In  this  work,  the  system  being  considered  is  a  pendulum,  where 
the  forcing  is  directed  along  a  line  at  some  angle  a  measured  with  respect  to  the  horizontal.  The 
dynamics  of  the  pendulum  are  given  by  the  second  order  o.d.e. 

16  -h  mtlv  cos{9  —  a)  +  mi,g£  sin  0  =  0,  (5) 

where  v{t)  =  r{t)  is  the  velocity  which  the  pendulum  hinge  point  is  forced  to  follow.  The  parameters 
are  nib  =mass  of  the  pendulum,  I  =the  planar  moment  of  inertia  of  the  pendulum  computed  about 
the  hinge  point,  and  t  =  the  distance  from  the  hinge  point  to  the  center  of  mass  of  the  pendulum. 
The  averaged  potential  for  the  mechanism  is 

^2/2 

Vyi(0)  = . cos^(0  -  a)  —  rntglcosOy  (6) 
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q  (t)  =pcosco^ 


Figure  1:  The  rapidly  forced  cart  and  n-penduluna  on  an  inclined  plane. 

where  =  (1/^)  lo  with  i;(‘)  being  a  piecewise  analytic  periodic  function  of  period  T. 

Although  the  exact  form  of  the  periodic  input  is  of  little  importance  (as  noted  in  [7]),  to  focus  the 
discussion,  suppose  i;(i)  =  Pusinujt^  corresponding  to  having  the  pendulum  hinge  point  execute  the 
motion  qo{t)  =  -jicosut  along  the  prescribed  axis.  In  the  nondimensional  version  of  this  problem 
(as  described  in  [18]),  two  dimensionless  parameters  play  a  role: 


.  _  _  9 

S=J  and  e  =  j^, 


in  terms  of  which  we  write  down  a  nondimensional  averaged  potential  Va{q)  as  described  in  the 
next  section.  In  rescaling  the  noncanonical  Hamiltonian  and  computing  the  averaged  potential,  it 
is  found  that  for  high  frequencies  (small  e),  the  averaged  system  tends  to  stabilize  in  the  direction 
of  oscillation.  Further  analysis  shows  that  the  nonautonomous  system  executes  periodic  motions, 
called  hovering  motions  [16,  18,  15],  about  the  averaged  fixed  point,  where  the  period  of  motion  is 
the  same  as  that  of  the  forcing,  and  the  amplitude  is  roughly  proportional  to  the  forcing  amplitude 
divided  by  the  forcing  frequency  squared  (=  Jc,  and  hence  the  amplitude  of  the  periodic  motions 
tends  to  zero  as  e  0).  Equilibria  of  the  averaged  system  and  periodic  orbits  of  the  nonautonomous 
system  are  shown  to  arise  as  the  result  of  bifurcations  at  critical  values  of  the  parameters  5  and  €. 
Furthermore,  it  has  been  shown  shown  that  there  is  a  close  correlation  between  the  averaged  phase 
space  and  Poincare  map  of  the  nonautonomous  system. 


In  this  report,  we  show  that  many  of  the  qualitative  features  of  the  single-degree-of-freedom 
(SDOF)  system’s  dynamics  also  appear  in  the  dynamics  of  periodically  forced  n-link  chains.  Ques¬ 
tions  concerning  charateristic  length  scales  become  more  important  and  more  subtle,  however. 
Specifically,  the  focus  of  this  report  is  the  bifurcation  and  stability  theory  of  a  periodically  forced 
cart  and  n-pendulum  on  an  inclined  plane,  illustrated  in  Figure  1.  Roughly  speaking,  the  question 
we  study  is  how  the  various  links  in  the  chain  align  themselves  when  the  cart  executes  high  frequency 
but  small  amplitude  periodic  motions.  As  with  the  SDOF  pendulum,  parts  of  this  story  have  roots 
dating  back  to  the  early  part  of  this  century [12].  Stabilization  of  the  inverted  n-pendulum  has  been 
more  recently  considered  by  Acheson[l].  In  this  work,  the  Floquet-theoretic  approach  is  applied  to 
a  linear  model  of  n-pendulum  dynamics  in  a  neighborhood  of  the  fully  inverted  equilibrium.  The 
main  result  of  this  report  is  a  bounding  of  the  forcing  amplitude  at  which  the  fully  inverted  equi- 
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librium  is  stabilized  given  the  minimum  and  maximum  natural  frequencies  {umin  and  u)rnax  in  [1]> 
resp.)  and  the  forcing  frequency  (cjo)-  Acheson  remarks  that  a  necessary  condition  for  stabilization 
is  cjn  >  Furhter  results  in  this  direction  will  be  discussed  below  in  Section  5. 

\J  TJitCLX 

An  advantage  of  taking  the  averaging  approach  outlined  in  this  introduction  is  that  it  provides 
the  means  to  find,  determine  the  stability,  and  study  bifurcations  of  all  equilibria.  While  it  is,  in 
principle,  possible  to  generally  follow  the  procedure  used  to  study  the  averaged  SDOF  system  as 
outlined  in  [18],  the  n-DOF  system  is  considerably  more  complex.  Specifically,  it  possesses  many 
more  symmetries  arising  from  mass  and  length  parameters,  whose  breakings  have  profound  effects 
on  the  equilibrium  and  bifurcation  structure  of  the  system.  The  phase  space  geometry  of  the  SDOF 
system  lends  itself  well  to  a  detailed  study  of  the  stability  and  basins  of  attraction  of  equilibria 
of  the  averaged  system  dynamics.  For  many-link  pendula,  the  picture  is  obscured  by  the  large 
number  of  dimensions  of  the  system’s  phase  space.  A  more  important  problem,  however,  is  that  in 
the  n-DOF  system  the  connection  between  the  averaged  system  dynamics  and  the  nonautonomous 
(controlled)  system  dynamics  depends  crucially  on  the  characteristic  length  scale  of  the  smallest 
system  component.  This  problem  is  computationally  tractable  for  the  double  pendulum,  whose 
dynamics  we  consider  in  Section  5.  In  the  general  case,  it  may  be  possible  to  use  the  averaged 
Hamiltonian  in  estimating  the  size  of  the  basins  of  attraction,  and  we  comment  on  this  possibility 
in  the  conclusion. 

2  The  Rapidly  Forced  Cart  and  n-Pendulum  on  an  Inclined 
Plane 

2.1  Averaging  and  Characteristic  Length  Scale  in  the  SDOF  Case 

Before  treating  the  general  case,  it  is  useful  to  describe  the  connection  between  characteristic  length 
scales  and  our  theory  of  averaging  in  the  case  of  a  SDOF  pendulum  on  a  cart.  The  goal  here  is  to 
understand  the  stability  characteristics  in  terms  of  the  critical  points  of  V>i(*)  given  in  (6).  These 
points  are  found  as  solutions  to  the  critical  point  equation 

21  g  sin  0  —  rribicr'^  sin  2{9  —  a)  =  0, 

and  when  a  =  |  (the  case  in  which  the  suspension  point  is  periodically  forced  along  the  vertical 
axis),  the  inverted  configuration  0  =  tt  is  always  a  critical  point.  The  averaging  theory  of  [4,  5]  states 
that  if  a  critical  point  is  a  local  minimum  of  V>i(*),  then  for  oscillatory  forcing  of  sufficiently  high 
frequency,  the  pendulum  will  execute  stable  motions  in  a  neighborhood  of  the  point.  In  the  present 
case,  the  second  derivative  V"(7r)  is  positive  precisely  when  >  gL  Making  the  apporoxima- 

tion  I  w  and  taking  the  forcing  function  v{t)  =  pujsin{ujt),  our  theory  predicts  that  there  will 

be  some  value  ojs  ~  such  that  the  pendulum  will  execute  stable  motions  very  close  to  the 

inverted  configuration,  provided  oj  >  This  recovers  the  classical  result  on  this  problem.  (C.f.  [2], 
p.  152.) 

Since  the  critical  frequency  parameter  value  of  the  averaged  potential  is  Ucr  —  it  is  of 

interest  to  know  how  closely  Ucr  approximates  cj^.  Since  iJcr  ~  ^  naive  conjecture  would 

be  that  the  stabilization  effects  of  oscillatory  forcing  are  observed  at  lower  frequencies  when  the 
pendulum  length  parameter  i  becomes  small.  Unfortunately,  this  conjecture  is  false,  and  we  have 
shown  in  [7]  that  Ug  ~  ojcr  only  when  ^3  ~  This  means  that  to  apply  averaged  potential  theory  at 
small  characteristic  length  scales,  we  need  to  have  ^  and  hence  Ug  ^  u)  cr 

For  the  n-pendulum  systems  we  treat  next,  there  are  two  characteristic  length  scales  which  are 
important — the  total  chain  length,  £,  and  the  length  of  the  shortest  link,  ig.  Of  the  two,  £s  is  the 
more  important  for  linking  the  critical  point  analysis  of  the  averaged  potential  to  the  stability  of 
the  forced  system.  Hence  in  the  rescaling  discussed  below,  this  length  scale  will  play  an  important 
role.  A  more  detailed  discussion  will  be  presented  in  Section  5. 
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2.2  System  Model 


Consider  the  system  illustrated  in  Figure  1,  where  each  rrii  is  a  point  mass  located  at  the  end  of  the 
z-th  massless  link,  li  is  the  length  of  the  z-th  link,  9i  is  the  angular  deviation  of  the  i-th  link  with 
the  hanging  vertical  position,  qo  is  the  position  of  the  cart  along  the  incline  of  angle  a,  and  g  is  the 
acceleration  due  to  gravity.  A  kinematic  analysis  of  this  system,  similar  to  that  presented  in  [14], 
results  in  the  derivation  of  a  Lagrangian  of  the  form  (1),  which  is  associated  with  the  noncanonical 
Hamiltonian 

n{e,P-,qo)  =  l{P-qoA{9)fMi9)-^iP-qoAm+Vaie-,qo),  (7) 

where 


Miji9) 

Ai{e;a) 

Va{e-,qo) 


y]  mjt  idj  cos{9i  -  9j) 
^fc=max(t,j)  J 

m*  j  cos{9i  -  a) 


U=1 


E 


j  gii{l  -  cosOi) 


t=l  L  \k=i 


(8) 

(9) 


We  immediately  note  that  H  is  dependent  on  the  variable  90  only  through  the  cart  velocity  go- 
Suppose  that  the  cart  position  go  is  given  explicitly  as  a  piecewise  analytic  periodic  function  of 
time — for  instance,  go(t)  =  /Scosut. 


2.3  Rescaling  and  Averaging 

Let  m  =  ‘rrii,  and  let  £s  be  the  smallest  link  length.  Then  the  length  of  all  other  links  can  be 
written  £i  =  Ai4,  where  1  <  Xi  <  k,  where  k  is  an  0(1)  constant.  We  nondimensionalize  (7)  by  the 
change  of  variables  qi  =  9i,  P  =  muitlp,  V.  =  and  a  change  of  time  scales  r  =  wt.  Let 

PQ{t)  =  go(t).  After  making  these  substitutions  and  simplifying,  we  obtain 


where 


H{q,p]T)  =  ^(p-^u(r)A(g))M(g)  ^(p-(5t;(T)A(g)) +V’o(g;T) 


MiM 

Ai{q]  a) 

Va{q;r) 

Vg{q) 


y  Pk  XiXj  cos(gi  -  qj) 

^fc=max(i,j)  / 

^  n  \ 

T  (E 


Kk=l 


E 

i=l 


—  cosgi) 


\k=i 


(10) 


(11) 


and  where  5  =  /3/4,  e  =  g/lgUj'^,  v{t)  =  ^Qir),  and  pi  =  rrii/ m.  The  averaged  Hamiltonian  H 
is  defined  as  H{q,p)  =  ^  H{q,p;T)dT,  where  T  is  the  fundamental  period  of  v{-).  The  averaged 
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Hamiltonian  has  the  general  form 

Hiq,p)  =  ^p^M{q)-^p  +  VA{q-,S,e)  (12) 

where  the  averaged  potential  VA{q\S^e)  is 

VAiq;S,e)  =  ^-^A{qfM{q)-^A{q)+eV,{q).  (13) 

Remark  1  The  parameters  S  and  e  have  very  clear  physical  interpretations.  5  represents  the  ratio 
of  the  forcing  magnitude  P  to  the  chain  length  t,  and  e  can  be  interpreted  as  the  ratio  of  the  system 
“natural  frequency”  Wn  =  squared  to  the  forcing  frequency  a;  squared.  In  assuming  small 

amplitude,  high-frequency  forcing,  we  are  making  the  assumption  that  P  and  Hence, 

0  <  5  1  and  0  <  e  1. 

Remark  2  An  advantage  of  using  the  nondimensional  model  is  that  it  facilitates  the  mathmematical 
discussion  concerning  equilibria  and  stability  in  several  physically  interesting  limits.  Of  particular 
interest  is  the  limit  in  which  the  forcing  frequency  becomes  infinite  while  the  forcing  amplitude 
remains  finite  but  non-zero.  In  the  nondimensional  model,  this  limit  is  taken  by  fixing  5  >  0 
while  letting  e  0.  Another  physical  interpretation  for  this  limit  is  that  the  usual  potential  forces 
(i.e.  gravity)  in  Va  become  negligible  when  compared  to  those  potential  forces  arising  from  the 
periodic  forcing.  While  the  theory  unfolds  most  easily  in  terms  of  the  e^5  parameters,  it  will  be 
useful  for  physical  interpretation  to  retain  also  a  dimensional  version  of  the  averaged  potential.  The 
dimensional  version  of  (13)  is 

V4(g;  <7)  =  ^A{efM{e)-^A{6)  +  V,{e\  (14) 

where 


with  T  being  the  fundamental  period  of  z;(-). 

Remark  3  For  the  sake  of  explicitness,  we  may  carry  out  the  nondimensionalization  of  the  chain 
model  assuming  sinusoidal  forcing  qo{t)  =  P  cosa;i,  but  the  assumption  of  this  particular  form  of  the 
periodic  input  is  not  necessary. 

3  Behavior  of  the  Averaged  System  in  the  Limit  e  — ^  0 

In  experiments,  it  is  usually  the  case  that  the  mechanism  producing  the  periodic  motion  of  the  pen¬ 
dulum  base  produces  motion  of  fixed  amplitude  and  variable  frequency.  This  constraint  results  in 
fixing  the  value  of  P  (or,  equivalently,  5),  while  allowing  us  to  vary  u)  (or  e).  Studies  of  stabilization 
for  SDOF  systems  have  shown  that  stabilization  is  a  high-frequency  phenomenon.  Therefore,  the 
approach  we  take  is  to  study  the  behavior  in  the  limit  as  e  ^  0  and  5  is  held  constant. 

We  begin  by  deriving  the  equilibrium  equations,  which  we  write 

M{q)^^p  =  0 

^l^^{A{qfM{q)-^A{q))+e^^  =  0.  (15) 

Stabilization  in  the  high-frequency  cj  oo  limit  is  easily  studied  by  simply  letting  e  =  0.  By 
definition,  is  nonsingular,  from  which  the  first  equation  of  (15)  implies  that  p  =  0  at  an 

equilibrium.  We  state  the  results  of  our  analysis  of  the  second  equation  in  the  following  theorems. 
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Theorem  1  Equilibria  associated  with  (10)  in  the  e  -4  0  limit  inelude  qi  =  a-\-  kTr/2,  where  kE.Ii. 


Proof:  Setting  e  =  0  in  (15)  leaves  us  with  n  equilibrium  equations  defined  by 

d 


dqi 


[A{qfM{q)-^A{q))=Q. 


(16) 


This  system  of  equations  may  be  written  in  symbolic  form  as 

\T 


where  is  a  diagonal  2-tensor  (matrix)  and  —  is  a  3-tensor.  Among  solutions  to  this 
system  are  those  for  which  one  of  the  following  conditions  hold: 


dq 


A{q)  =  0 
dM{q)^^ 


dq 


=  0, 


(17) 

(18) 


since  nonzero  components  of  of  the  form  sin(gi  -  a),  while  nonzero  components  of  —  are 
of  the  form  sm{qi  -  qj).  This  proves  the  theorem.  □ 


Theorem  2  Equilibria  q^  of  (15)  with  coordinate  entries  qi  =  a  +  1  <  i  <  n  and  ki  £  Z 

are  stable  rest  points  of  the  averaged  system. 

Proof:  We  prove  the  theorem  by  showing  the  Hessian  matrix  H(g)  =  ^  Qg2^^  is  positive-definite 
for  each  qe  of  the  prescribed  form.  Positive-definiteness  implies  that  qe  is  a  local  minimum  of  the 
averaged  potential  Va(^),  and  therefore  qe  is  a  stable  rest  point  of  the  averaged  Hamiltonian  system. 
Computing  partial  derivatives,  we  obtain 


H(5) 


+ 


.d^Ajqy 

dq^ 


■M(g)-^  +3 


^dA{q)'^  dM{qy^ 


dq 


dq 


+  24(g)' 


dq^  , 


At  qs,  recall  that  A{qs)  =  0,  and  therefore  (19)  simplifies  to 

For  any  vector  x  E  K” ,  we  have 


q-Qs 


24(g) 


(19) 


x  =  y^M{qs)  ^y>0, 


q-Qs 


where  we  have  substituted  y  =  x.  Because  M(g)'"^  is  positive  definite  by  definition, 

H(gs)  is  also  positive  definite.  □ 


Theorem  3  Equilibria  q^  =  a  +  kiTT  experience  a  bifurcation  at  e  =  0.  Further,  there  exists  ei  -C  1 
such  that  for  0  <  e  <  ci,  these  equilibria  are  unstable. 
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Stable 

Equilibria 


Unstable 

Equilibria 


Figure  2:  Selected  equilibria  for  the  averaged  cart  and  double  pendulum  in  the  limit  as  e  0.  As 
stated  in  Theorem  1,  equilibria  aligned  with  the  direction  of  motion  are  stable,  and  those  aligned 
normal  to  the  direction  of  motion  are  unstable. 


Proof:  At  gi  =  a  -1-  nii:,  both  =  0  and  -  =  0,  in  which  case  H(g)  simplifies  to 


dq 

-1 


=  )-■..(,)-«  I  AW. 


(20) 


We  can  eliminate  either  or  by  computing  another  partial  derivative  of  (18):  solving 

we  obtain 

02M(g)-^ 


for  ^  -  and  setting  ^4^  =  0  and  —  =  0,  we  obtain 


Qgi  •  and  setting  —  u  anu  — ^ 

d^A{q)  _  _1 
2 


=  17  M{q). 


dq^  ~  dq^ 

We  substitute  this  last  expression  into  (20)  and  find 

5g2 


H(g)  =  (  ]  A{q)  =  0. 


5^2 


-) 


□ 


Remark  4  Theorem  1  indicates  that  in  the  e  0  limit,  the  system  stabilizes  in  configurations  where 
the  links  are  aligned  with  the  cart  motion.  There  are  unstable  equilibria  prescribed  by  configurations 
where  the  links  are  aligned  with  the  direction  normal  to  the  cart  motion.  This  is  illustrated  for  the 
cart  and  double  pendulum  in  Figure  2.  We  consider  the  equilibria  and  bifurcations  of  the  cart  and 
double  pendulum  in  Section  5. 


4  Bifurcations,  Stability,  and  Behavior  of  the  Vertically  Forced 
n-Pendulum  as  n  — >  oo 

In  this  section,  we  consider  stabilization  as  we  pass  to  the  limit  in  which  we  let  n  ->  oo  while  holding 
the  pendulum  length  and  mass  constant.  The  main  result  to  emerge  from  this  section  is  that  as  we 
add  degrees  of  freedom,  the  oscillation  frequency  at  which  stabilization  occurs  becomes  large,  and 
ultimately  tends  to  infinity. 


4.1  Bifurcations  of  the  Vertically  Forced  n-Pendulum 

For  the  vertically  forced  pendulum  systems  (a  =  |),  the  link  configurations  in  which  all  qi  are  either 
0  or  TT  are  stationary  points  (equilibria)  of  the  forced  Hamiltonian  system  (7)  for  all  values  of  the 
input  qo{')^  (Cf.  Theorem  1  and  Theorem  2.)  Although  these  points  are  always  equilibria,  their 
stability  characteristics  change  as  we  vary  the  frequency  of  the  oscillatory  forcing.  This  is  reflected 
in  bifurcations  of  the  averaged  potential  described  as  follows. 

Theorem  4  In  the  coefficient  system  (11)  associated  with  the  nondimensionalized  Hamiltonian, 
suppose  a  =  |.  Let  denote  the  set  of  critical  points  of  the  averaged  potential  VA{q)  defined 
in  (13)  in  which  each  qi  is  either  0  or  n,  1  <  i  <  n.  For  all  qe  G  {qs}  except  the  trivial  hanging 
equilibrium  =  (Oj » •  •  ?  0)^?  there  exists  a  critical  value  of  e  denoted  Ccr  >  0  such  that  at  e  =  Ccr,  qe 
experiences  a  bifurcation  such  that  for  0  <  e  <  Ccr  7  qe  a  strict  local  minimum  of  Va- 

Remark  5  Recall  that  e  ^  where  u  is  the  frequency  of  the  oscillatory  forcing.  Theorem  4  states 
that  for  each  equilibrium  in  which  all  links  are  aligned  vertically,  there  is  a  critical  frequency  u)cr 
above  which  the  equilibrium  is  a  relative  minimum  of  the  averaged  potential 


Proof  of  Thm.  4:  Prom  Theorem  2,  equilibria  belonging  to  the  set  are  stable  for  the 
averaged  system  when  e  =  0,  so  what  remains  to  be  done  to  prove  the  theorem  is  to  show  that 
for  each  qe  G  {qs}  there  exists  an  Ccr  >  0  at  which  a  bifurcation  occurs,  changing  the  stability 
type  of  the  critical  point.  Recall  that  the  stability  of  a  critical  point  of  (13)  may  be  determined 
from  the  concavity  of  the  averaged  potential  ^^(g),  which  is  inferred  from  the  Hessian  matrix  H(g). 
Furthermore,  the  equilibrium  experiences  a  bifurcation  if  for  some  set  of  parameters  the  averaged 
potential  becomes  loses  rank,  i.e.  if  the  Hessian  matrix  has  one  or  more  zero  eigenvalues.  For 
9e  ^  {gs}j  H(g)  can  be  written 


H(9e)  =  ( 


(52  dA{qy 

2  dq 


■Miq)- 


.,dA{q)  d^V,{q) 


dq 


—  e- 


9g2 


?=9. 


Recall  from  the  proof  of  Theorem  2  that  the  matrix  is  always  positive-definite, 

and  furthermore  is  symmetric  (because  M{q)~^  is  symmetric  and  positive  definite).  By  definition  of 
^s(5')i  is  a  diagonal  matrix  which  is  easily  partitioned  into  a  positive-definite  part  we  shall 

denote  V(5)+,  and  a  negative-definite  part  denoted  -V~(q)  where 'V(q)~  is  itself  a  positive-definite 
matrix.  Note  that  for  the  trivial  hanging  equilibrium  that  V“(g)  =  {0}"^”,  in  which  case  the 
Hessian  is  always  positive-definite  and  the  equilibrium  is  always  stable.  For  all  other  equilibria,  the 
Hessian  matrix  can  be  written  as  the  difference  of  positive-definite  matrices  as  follows: 

-‘[VWIU.- 

9=9. 

The  theorem  will  be  proved  by  establishing  the  following  lemma. 


H(9e)  = 


(52  dA(q)^ 


dq 


M{q) 


-idAjq) 

dq 


+  eY{qy 


Lemma  1  Let  Pi  and  P2  be  n  x  n  positive-definite  symmetric  matrices.  Then  there  exists  a  real 
number  €c  >  0  such  that  Pi  —  eP2  is  positive  definite  for  all  e  <  Cc  and  Pi  —  eP2  is  not  positive 
definite  for  e>  Cc* 


Proof  of  lemma:  Pi  “  r}P2  is  symmetric  for  all  values  of  the  parameter  77,  and  because  Pi  is 
positive  definite,  for  any  vector  a:  ^  0  of  the  appropriate  dimension,  x'^Pix  >  0.  Given  such  an  x, 
it  is  clear  that  there  is  some  number  771  >  0  such  that  x’^Pix  —  7}x'^P2^  <  0  for  all  77  >  771.  Indeed, 
given  X,  we  can  chose  771  =  771  (x)  =  x^Pix  /x^P2X.  The  that  we  seek  is  the  inf  of  771  (x)  over  all 
X  such  that  ||x||  =  1.  O 
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Making  the  appropriate  substitutions  for  Pi  and  P2  in  the  previous  lemma  and  letting  Ecr  be  the 
smallest  €c>  there  exists  for  each  G  {Qs}  an  Ccr  >  0  at  which  H(ge)  is  no  longer  positive-definite. 
Given  that  for  qe  G  {qs},  the  terms  ^  M {q)~'^  ,  V(g)'^,  and  V(g)“  are  constant  with  re¬ 

spect  to  q  and  that  the  term  eV(g)'”  grows  linearly  with  e,  we  can  conclude  that  H(ge)  is  never 
positive  definite  for  e  >  Ccr-  Hence,  each  q^  is  rendered  unstable  by  a  bifurcation  at  some  e  =  Ccr  >  0 
as  e  increases  in  value.  □ 


Remark  6  Note  that  the  type  of  bifurcation  experienced  by  equilibria  in  the  set  {qs}  must  necessar¬ 
ily  be  pitchfork,  in  that  exactly  two  equilibria  must  bifurcate  symmetrically  about  elements  of  {qs}- 
Note,  however,  that  this  does  not  imply  that  other  types  of  bifurcations  cannot  occur  elsewhere. 
In  fact,  in  our  study  of  the  double  pendulum  in  the  next  section,  we  shall  see  that  saddle-node 
bifurcations  do  occur  away  from  elements  of  {g^}. 


4,2  Stabilization  of  the  Fully  Inverted  Equilibrium  as  n  — 00 

Theorem  5  Consider  the  fully  inverted  equilibrium  gg  =  (tt,  . . .  ,7r)^  of  the  n-linkj  vertically  {a  = 
7r/2)  forced  chain.  Assume  that  all  masses  and  link  lengths  are  equal,  and  that  the  total  mass  and 
length  of  the  chain  is  held  constant;  i.e.  fii  =  ~  and  Xi  =  1  for  1  <  i  <  n.  Let  Hn(ge)  denote 
the  Hessian  matrix  associated  with  this  equilibrium,  and  let  Ccr  be  the  smallest  value  of  e  for  which 
detHn(ge)  =  0.  Then  Ccr  is  strictly  decreasing  as  n  00,  and  the  critical  frequency  ujcr  o.t  which 
the  equilibrium  is  stabilized  in  the  averaged  system  dynamics  tends  to  infinity  (00)  as  n  infty. 


Proof;  Our  strategy  of  proof  is  to  convert  the  determinant  calculation  to  an  eigenvalue  problem, 
where  the  eigenvalues  we  compute  correspond  to  critical  values  of  e  and  the  smallest  of  these  values 
is  Ecr-  Because  of  the  special  structure  of  Hn(ge),  we  will  be  able  to  show  that  Ecr  for  the  n-DOF 
system  is  strictly  smaller  than  for  the  (n  —  l)-DOF  system  for  all  n.  We  begin  by  explicitly  writing 
out  the  determinant  as 


det  Hn(ge)  =  det 


dAjq)'^ 
2  dq 


M{q) 


-idAjq)  ,  ^d^Vgjq) 


dq 


-h  6- 


9g2 


(21) 


We  now  proceed  to  write  each  of  the  terms  in  the  preceding  equations  in  terms  of  n.  First,  we  note 
that  SAidl  is  a  diagonal  matrix  of  the  form 


dA{q) 

dq 


0  -(Er=2Mi)^2 


0 


0 


/  -n 


0 


0  \ 


0  -(n  — 1) 


VO  0  -1  / 

Similarly,  --^1-  can  be  reduced  to  a  diagonal  matrix  of  the  form 


dq^  n 


(  -n  0 
0  -(n  -  1) 


V  0 


0  \ 
0 

-ly 


0  \ 
0 

f^n^n  ) 


-A. 

n 


=  iv. 

n 
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For  we  first  consider  a  reduction  of  M{q):  specifically,  we  have 


M{q)  = 


1 


n 


/  n  n-l  n-2 

n~l  n-l  n-2 

n  —  2  n  —  2  n  —  2 


V 


1  1 

The  inverse  of  M  given  in  (22)  can  be  expressed  as 


M-i  = 


1\ 
1 
1 

ly 


=  -M. 

n 


(22) 


1 

-1 

0  ••• 

0 

0 

0 

-1 

2 

_1  ... 

0 

0 

0 

0 

-1 

2  ••• 

0 

0 

0 

0 

0 

0 

2 

-1 

0 

0 

0 

0  ••• 

-1 

2 

-1 

0 

0 

0  ••• 

0 

-1 

2 

y 

for  all  n.  Then  M{q)~^  —  nM“^,  and  then  (21)  can  be  written 

detH„(5e)  =  deti  +  =  det-Hn(ge)  =  0, 

n  \  2  /  n 


(23) 


where 


j  f»(^n—  «)  — ^n(n  —  1)  0 

-^n{n-l)  2(n  -  1)  (^(n  -  I)  -  <)  -  ^  (r.  -  l)(n  -  2) 

0  2)  2(n  -  2)  (^(n  -  2)  -  €) 


V 


-4 


1  (^2  -  c) 

-4  2(4-«)  ) 


H„(ge)  is  a  Jacobi  matrix.  The  theorem  can  be  proven  by  converting  (23)  to  an  eigenvalue  problem 
and  making  use  of  the  following  lemma  for  Jacobi  matrices,  which  is  given  in  [11]. 


Lemma  2  Let  L„  be  ann  x  n  Jacobi  matrix  of  the  form 


L„  = 


/  6i  Cl  0 

02  i>2  C2  0 

0  03  63  C3 


0  \ 
0 
0 


0 
V  0 


^n— 1  1  ^n— 1 

0  ttn  J 


where  ai,  hi^a  G  E  for  all  i  and  let  Lr  denote  the  principal  r  xr  submatrix  Ln[l, . .  -,^11)  •  •  •  j^]-  U 
aiCi-i  >  0,  i  G  2, . . .  ,n,  then 

1.  all  characteristic  roots  o/Ln  are  real  and  simple,  and 

2.  between  any  two  characteristic  roots  of  Ln  lies  exactly  one  characteristic  root  of  Ijn-i  • 


Proof:  This  lemma  appears  in  [11],  and  its  proof  may  be  found  in  [3]. 


O 
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We  now  premultiply  Hn(?e)  by  a  positive-determinant,  n  xn  diagonal  matrix  C  and  postmultiply 
by  a  positive-determinant,  nxn  diagonal  matrix  Q,  where 


1 

Tl 

0 

0 

f  1 

0 

0  ••• 

0  \ 

0 

1 

n— 1 

0  ••• 

0 

0 

1 

2 

0 

0 

0 

0 

1 

n-2 

0 

,  Q  = 

0 

0 

1 

2 

0 

(24) 

0 

0 

0  ••• 

1 J 

u 

0 

0  ••• 

1 ) 

to  obtain  _  _ 

det  [CH„(ge)Q]  =  detH„(ge)  =  0, 


where 


1 

1  - 

O 

0 

0  \ 

^(n-l)-e 

-f(n-2) 

0 

0 

0 

-f(n-l) 

f  (n-2)  -€ 

0 

0 

0 

0 

0 

•••  <52 -e 

4 

0 

0 

0 

2 

Equation  (25)  is  now  in  the  form  of  the  eigenvalue  problem 


(25) 


detHn(ge)  =  det(Ln  -  ^I)  =  0, 


where  is  a  positive-definite  Jacobi  matrix  of  the  form  considered  in  Lemma  2.  Note  that  because 
Ln  is  positive-definite,  min(e)  >  0.  We  now  note  that  we  can  obtain  the  analogous  matrix  L^-i  for 
an  n  —  1-DOF  system  by  blocking  off  the  first  row  and  column  of  Ln.  Denote  the  eigenvalues  of 
Ln.-i  as  {cn-i}  and  the  eigenvalues  of  Ln  as  {cn}*  By  Lemma  2,  the  eigenvalues  of  Ln-i  lie  between 
those  of  Ln.  Hence,  we  have 

0  <  min({en})  <  min({en-i}) 

for  all  n.  Therefore,  as  n  oo,  min({en})  -)•  c,  c  >  0.  To  complete  the  proof  of  the  theorem,  we 
note  that  the  critical  frequency  is  related  to  e„  by 


€n  = 


Since  the  length  of  the  shortest  link  is  ~  1/^j  the  result  of  the  theorem  follows. 


□ 


Remark  7  Theorem  5  implies  the  critical  frequency  Ucr  at  which  each  element  of  {^s}  is  rendered 
stable  increases  as  additional  links  are  added  to  the  pendulum  (while  holding  the  total  mass  and 
length  of  the  pendulum  constant).  Consequently,  in  the  limit  as  the  number  of  links  n  — oo, 
Ucr  — ^  00  for  each  element  of  {qs}^  This  conclusion  suggests  that  it  is  therefore  impossible  to 
stabilize  any  equilibria  in  the  set  in  the  continuum  limit;  i.e.  it  is  impossible  to  stabilize  an 
inverted  (planar)  string. 

Theorem  5  can  be  easily  verified  using  widely  available  computer  algebra  software.  In  this  study, 
we  compute  the  eigenvalues  of  the  nxn  matrix  L^  for  1  <  n  <  50  and  6  =  0.2.  The  results  of  this 
calculation  are  shown  in  Figure  3,  where  we  see  in  the  left  frame  that  min(^n)  tends  toward  zero  as 
n  becomes  large.  Consequently,  0  quickly,  as  shown  in  the  right  frame. 
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n  n 

Figure  3:  det(Ln)  vs.  n  (left)  and  ter  vs.  n  for  the  fully  inverted  equilibrium  of  the  vertically  forced 
chain.  As  predicted  by  Theorem  5,  both  det(Ln)  and  Ccr  decrease  as  n  increases. 

4.3  Stabilization  of  the  Completely  Folded  Equilibrium  as  n  oo 

By  employing  essentially  the  same  procedure  used  to  prove  Theorem  5,  we  can  prove  that  in  the  limit 
as  n  00,  Ccr  tends  to  zero  for  any  regular  pattern  of  link  configurations.  In  the  next  section,  we 
present  numerics  which  suggest  that  the  fully  folded  configuration  is  one  of  the  first  configurations 
to  stabilize  as  e  0,  so  it  is  reasonable  to  expect  that  values  of  Ccr  for  those  configurations  also 
tend  to  zero. 

Remark  8  In  [18]  it  was  seen  in  the  SDOF  system  that  for  u  <  Uen  the  averaged  phase  portrait 
generally  reflected  the  Poincare  map  of  the  nonautonomous  system,  specifically  in  terms  of  peri¬ 
odic  orbits  associated  with  averaged  fixed  points  and  the  stability.  Periodic  orbits  associated  with 
averaged  fixed  points  manifest  themselves  as  1:1  resonances  in  the  Poincare  map.  What  the  aver¬ 
aged  phase  portrait  never  captures  are  subharmonic  resonant  responses  in  the  Poincare  map.  For 
uj  >  Ucr,  it  was  observed  in  [18]  that  these  subharmonic  resonances  do  not  play  a  significant  role  in 
the  nonautonomous  dynamics.  This  has  also  been  established  rigorously  in  two  limiting  cases  of  the 
parameters  S  and  e.  For  uj  <  ujcr,  however,  subharmonic  responses  do  play  a  significant  role  in  the 
nonautonomous  dynamics,  along  with  homoclinic  tangling  of  phase  space  separatrices  and  the  asso¬ 
ciated  bands  of  stochasticity.  It  is  therefore  likely  that  subharmonic  resonances  also  play  a  significant 
role  in  the  dynamics  of  the  periodically  forced  n-pendulum  for  uj  <  cjcr,  and  in  the  “large  n”  limit 
dominate  the  dynamics  for  all  u.  For  physical  systems,  damping  and  other  dissipative  forces  (e.g. 
aerodynamic  drag)  tend  to  mitigate  the  chaotic  behaviors  exhibited  by  purely  conservative  systems, 
and  one  would  expect  physical  systems  to  exhibit  fairly  regular  periodic  response  to  forcing. 

5  Numerical  Studies 

5.1  Bifurcations  of  the  Cart  and  Double  Pendulum 

The  results  from  Section  2  capture  the  equilibrium  behaviors  of  the  averaged  cart  and  n-pendulum 
in  the  case  of  infinite  forcing  frequency,  which  we  express  in  our  nondimensional  model  by  allowing 
e  0.  As  e  is  perturbed  off*  zero,  the  equilibria  found  in  Theorem  1  move  as  the  contribution 
of  gravity  to  potential  energy  becomes  significant,  and  for  distinct  critical  values  of  e  the  system 
experiences  bifurcations.  In  this  section,  we  focus  on  obtaining  numerical  results  for  the  vertically 
forced  double  pendulum  and  comment  on  bifurcations  for  the  symmetric  {fii  =  /i2  =  Ai  =  A2  =  1/2) 
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Figure  4:  Bifurcation  diagram  for  the  averaged  cart  and  double  pendulum,  where  /^i  =  /X2  =  Ai  = 
A2  =  1/2,  a  =  7r/2,  and  6  =  0.2.  In  the  plots,  we  see  that  the  inverted  and  folded  equilibria  are 
rendered  stable  in  transcritical  pitchfork  bifurcations.  These  bifurcations  give  rise  to  the  unstable 
equilibria  which  are  normal  to  the  direction  of  motion. 


system. 

The  necessary  equilibrium  equations  and  Hessian  matrix  are  easily  derived  from  the  general  form 
of  the  averaged  potential  given  in  (13).  The  reader  is  referred  to  [13]  for  full  details  of  this  calcula¬ 
tion.  We  now  turn  our  attention  to  calculations  of  equilibria  and  stability  for  the  vertically  forced, 
symmetric  (^1  =  =  Ai  =  A2  =  1/2)  double  pendulum,  where  we  fix  5  =  0.2.  The  algorithm 

which  we  use  is  a  two-dimensional  bisection  algorithm  which  utilizes  the  fact  that  the  equilibrium 
equations  for  this  system  define  two-dimensional  hypersurfaces  on  the  torus.  Equilibria  are  then 
intersections  of  the  intersections  of  each  surface  with  zero. 

For  the  vertically  forced  (a  =  7r/2)  case,  we  obtain  the  bifurcation  diagram  shown  in  Figure  4.  In 
this  figure,  we  see  that  for  large  values  of  e  the  system  possesses  unstable  equilibria  in  configurations 
aligned  with  the  vertical.  As  e  0,  new  equilibria  are  created  in  a  series  of  transcritical  pitchfork 
bifurcations  off  the  vertically  aligned  equilibria.  These  bifurcations  have  the  effect  of  changing  the 
stability  of  the  inverted  and  folded  equilibria  from  unstable  to  stable,  and  in  the  process  create  the 
unstable  equilibria  which  in  the  e  0  limit  correspond  to  the  configurations  normal  to  the  direction 
of  motion. 


5.2  Stabilization  Sequence  for  Vertically  Forced  Chains 

As  the  frequency  parameter  e  decreases,  the  averaged  potential  for  the  n-link  chain  undegoes  a  rich 
set  of  bifurcations,  and  many  equilibria  are  stabilized  by  the  oscillatory  forcing.  Some  indication  of 
this  richness  is  provided  in  Figure  5  where  for  two,  three,  and  four  link  chains  we  give  the  stabilization 
sequence  of  equilibria  in  which  successive  links  have  the  displayed  “up- down”  patterns.  Here  the 
bifurcation  parameter  is  the  r.m.s.  value  of  the  forcing  a  =  When  cr  =  0,  only  the  hanging 

equilibrium  is  stable,  and  this  is  the  only  local  minimum  of  the  averaged  potential.  As  a  increases 
(corresponding  to  an  increase  in  amplitude  or  frequency  of  the  forcing),  however,  the  form  of  the 
averaged  potential  changes  and  other  critical  points  become  local  minima.  There  are  a  number  of 
interesting  features  in  these  patterns,  including  the  fact  that  the  case  in  which  all  links  are  inverted 
is  the  last  to  become  a  relative  minimum  of  the  averaged  potential  (as  the  forcing  frequency  is 
increased)  for  chains  with  any  number  of  links. 
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Figure  5:  Bifurcation  Sequence  of  Heavy  Chains 


Another  observation  worth  noting  is  how  the  bifurcations  of  critical  points  of  the  averaged  po¬ 
tential  depend  on  other  physical  parameters.  Figure  6  illustrates  how  the  number  of  stable  equilibria 
depends  jointly  on  characteristic  link  lengths  together  with  r.m.s.  value  a.  In  this  figure,  the  various 
curves  define  boundaries  in  the  (^,  parameter  space  at  which  bifurcations  occur.  Only  the  first 
quadrant  of  the  (^,a^)-plane  is  physically  menaingful,  and  this  is  partitioned  into  four  sectors.  In 
each  sector,  it  is  precisely  those  configurations  depicted  by  the  cartoons  whose  stabilization  (under 
the  oscillatory  forcing)  is  predicted  by  the  analysis  of  the  averaged  potential.  Note  that  it  is  only 
in  the  topmost  sector  that  stability  of  all  vertically  aligned  configurations  is  predicted.  Here  what 
we  see  is  that  the  links  with  shorter  lengths  tend  to  become  stabilized  in  the  vertically  upright 
configuration  before  the  links  with  longer  lengths.  Clearly,  if  we  consider  chains  with  nonuniform 
links,  the  bifurcation  patterns  become  very  rich  as  the  number  of  links  increases. 

These  figures  raise  the  question  of  how  precisely  the  averaged  potential  analysis  predicts  the 
stabilization  of  an  n-pendulum  under  oscillatory  forcing.  To  some  extent  this  remains  an  open 
question.  Consider  the  dimensional  version  of  the  averaged  potential 

v.4(g)  =  +  V(g). 

As  u(-)  ranges  over  a  set  of  continuous  periodic  functions  of  period  1  with 

11^112=  [  v{tfdt  =  0^, 

Jo 

let  =  r}v{rit).  Then  the  r.m.s.  value  of  is 

Jo 

and  for  a  fixed  value  of  /?,  (as  u(-)  ranges  over  a  family  of  periodic  functions  of  L2-norm  /3,  we  ask 
how  much  the  stabilizing  effects  of  oscillatory  forcing  by  Vtj{‘)  can  vary. 

This  question  has  been  answered  explicitly  for  a  one-paerameter  family  of  piecewise  continuous 
inputs  in  the  case  of  a  single  DOF  pendulum  in  [7].  It  is  useful  to  summarize  this  result  before 
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Figure  6:  Bifurcation  locus  of  2-link  nonuniform  chains.  The  various  curves  define  boundaries  in  the 
(^,  (7^)-parameter  space  at  which  bifurcations  occur. 


addressing  the  n-DOF  case.  Take  u  to  be  a  rectangular  wave: 


i) 


h 


if  0  <  t  <  It; 
w  <t  <  1 


whose  definition  we  extend  to  all  t  G  [0,  oo)  by  time-1  periodicity:  v{t  H- 1)  =  v{t).  The  condition 
that  V  is  also  periodic  is  satisfied  provided  that  k  =  wh/{l  —  w).  The  dynamics  of  the  system  and 
the  averaged  potential  are  given  by  (5)  and  (6)  respectively  with  oc  =  To  study  the  dependency 
on  frequency  of  the  forcing,  let  Vrj{i)  =  T]v{r]t).  The  corresponding  value  of  is 


a 


2 

V 


12 


The  averaged  potential  analysis  of  the  stability  of  the  “inverted”  equilibrium  proceeds  as  follows. 
For  all  7],6  =  7r  satisfies  the  critical  point  equation  V^(0)  =  0,  and  since  V^(7r)  =  / (121)- 

rrihgi,  tt  is  a  strict  local  minimum  of  Va  whenever  rj  >  ^/12Ig/{mbihw).  Consider  the  linearized 
dynamics  of  (5)  in  a  neighborhood  of  {0,6)  =  (7r,0): 


16  —  TTibivO  —  rribgiO  =  0. 


(26) 


(See  [5]  or  [15]  for  information  on  linearized  models  of  Lagrangian  control  systems.)  We  shall  apply 
a  Floquet  type  analysis  to  the  dynamics  of  (26)  in  terms  of  the  time-T  mapping.  (T  =  I/77  is  the 
period  of  the  forcing.)  Specifically,  letting  {xk,yk)  =  i6{k/r]),6{k/r})),  this  map  is  given  explicitly 
by 

/  \  /  \ 

(27) 


where 


Xk+l 

yk-\-i 

1  —  w 


ill) 


0 
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and 


_w(  0  1 \ 

^(g  +  ri^h)  0  )■ 

Normalizing  the  rectangular  wave  so  that  hw  =  1,  the  second  derivative  of  the  averaged  potential 
becomes  positive  when  r]  >  rjcr  =  y/l2Igl{TnbC).  The  averaging  theory  of  mechanical  systems 
([5])  asserts  that  there  is  some  value  tj*  such  that  for  all  r)  >  rjs  the  linearized  model  (26)  will  be 
Lyapunov  stable.  Within  the  class  of  rectangular  wave  forms  we  have  specified,  rjcr  will  in  fact  vary 
even  though  the  r.m.s.  value  a  on  which  the  averaged  potential  depends  remains  constant.  More 
specifically,  when  hw  =  1,  we  shall  hold  =  ‘q/12  constant.  Then,  letting  w  range  between  0  and 
1,  the  shape  of  the  forcing  function  u(-)  changes  from  a  large-amplitude,  short-duration  up-pulse 
followed  by  a  relatively  shallow  negative  value  (when  ru  ~  0)  to  a  relatively  long  shallow  positive 
value  followed  by  a  large-amplitude  down-pulse  (when  u;  ~  1).  This  difference  in  wave-form  shape 
influences  the  range  of  values  of  vy  for  which  (26)  is  Lyapunov  stable.  The  analysis  is  carried  out 
in  terms  of  the  time-T  map,  (27),  and  in  particular,  we  conclude  that  (26)  is  stable  in  the  sense  of 
Lyapunov  precisely  when  the  eigenvalues  of  lie  on  the  unit  circle.  (For  more  details  on  this 

Floquet  type  of  stability  analysis,  see  [5].) 

For  a  single  DOF  pendulum,  it  is  not  difficult  to  find  a  fairly  explicit  relationship  between  and 
pulse- width  parameter  w.  In  [7]  we  make  the  somewhat  surprising  observation  that  ry*  does  not  vary 
a  great  deal  as  a  function  of  w,  although  a  waveform  characterized  by  a  prominent,  large-amplitude 
downward  pulse,  with  relatively  slow  return  upwards  (i.e.  in  terms  of  the  parameter  w  we  have 
w  ~  1)  will  stabilize  the  system  at  a  lower  frequency  than  waveforms  for  which  w  is  not  close  to  1. 
The  main  result  of  [7]  is  that  in  comparing  the  bifurcation  value  rjcr  with  the  critical  stabilization 
frequency  rjs  we  have  ry*  ~  T]cr  provided  ~ 

For  a  more  detailed  picture,  consider  a  pendulum  with  parameter  values  mb  =  1/2,  I  =  mbi^, 
and  take  g  =  10.  To  examine  the  case  in  which  0  is  scaled  in  direct  proportion  to  the  charateristic 
link  length  we  take  P  =  £.  Then  rjcr  =  ^120/^,  and  in  Figure  7  we  plot  rj  =  rjs /rjcr  as  a  function 
of  the  pulse-width  parameter  w.  We  again  see  that  the  stabilizing  effect  of  the  oscillatory  forcing 
is  relatively  insensitive  to  variation  in  the  waveform.  The  more  remarkable  feature  of  Figure  7, 
however,  is  that  the  relationship  between  rjs  and  rjcr  is  totally  independent  of  the  characteristic 
length  scale  £. 

The  stabilization  effects  of  oscillatory  forcing  of  n-DOF  pendula  depend  more  subtly  on  charac¬ 
teristic  length  scales,  and  it  can  be  shown  that  the  forcing  amplitude  should  scaled  with  the  shortest 
link  length  (corresponding  to  the  fastest  time  constant  in  the  system)  in  order  to  ensure  that  the 
key  qualitative  features  of  the  system  response  are  predicted  by  analyzing  the  averaged  potential. 
To  understand  this  remark  in  a  bit  more  detail,  we  consider  a  simple  two-link  system.  If  we  vary  the 
relative  link  lengths  in  the  chain  (as  depicted  by  the  cartoon  systems  in  Figure  6),  the  relationship 
between  the  bifurcation  values  rjcr  of  the  averaged  potential  and  the  actual  critical  stabilization 
frequencies  can  change  a  great  deal.  This  is  illustrated  in  Figure  8  where  we  plot  rj  —  rjs /rjcr  as  a 
function  of  the  pulse-width  parameter  w  for  the  configuration  (^1,^2)  =  (^TjO)  in  the  case  of  link 
lengths  £i  =  0.9  and  ^2  =  1-  Here  it  is  no  longer  the  case  that  the  graph  remains  unchanged  as  £i 
varies,  and  in  fact  we  have  found  that  the  slope  of  rj{w)  becomes  steeper  as  £i  is  decreased  rela¬ 
tive  to  £2  =  1.  For  this  simple  family  of  rectangular  wave-form  inputs,  the  stabilizing  effect  of  the 
periodic  forcing  is  increasingly  sensitive  to  the  pulse-width  parameter  as  the  relative  link  lengths 
become  more  disproportionate.  This  suggests  there  may  be  a  governing  robustness  principle  wherein 
increasing  the  degree  of  nonuniformity  in  the  pendulum  system  will  make  it  more  difficult  to  apply 
averaged  potential  theory  to  determine  stability  conditions  for  the  forced  system. 

5.3  Simple  Experiments  in  Open-Loop  Control 

Experiments  in  open-loop  oscillatory  control  for  pendulum  systems  are  not  terribly  difficult  to  con¬ 
struct  and  yield  compelling  demonstrations  of  the  effects  of  high-frequency  forcing  in  the  stabilization 
of  mechanical  systems.  The  apparatus  with  which  we  performed  our  experiments  was  constructed 
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Figure  7:  For  the  single  DOF  (simple)  inverted  pendulum,  the  normalized  stabilization  frequency 
7]  =  Vs/rjcr  is  plotted  (solid  line)  as  a  function  of  pulsewidth  parameter  w.  The  graph  remains 
unchanged  for  all  values  of  pendulum  length  £. 


Figure  8:  For  the  2-DOF  pendulum,  this  graph  depicts  the  normalized  stabilization  frequency 
rj  =  ris  frier  for  the  configuration  {61^62)  =  (tt,  0).  77  is  plotted  (solid  line)  as  a  function  of  pulsewidth 
parameter  w.  The  link-length  parameters  are  £i  =  0.9  and  £2  =  1.  This  graph  changes  with  the 
average  slope  of  the  solid  line  becoming  more  steep  as  £i  is  reduced  relative  to  the  length  £2  =  1. 
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Figure  9:  View  of  the  mechanism  which  rectifies  rotary  motion  to  periodic  vertical  motion. 


from  readily  available  plexiglass  and  plywood  stock,  with  a  common  toy  d.c.  motor  and  strips  of 
aluminum.  The  mechanism  which  rectifies  the  rotary  motion  of  the  motor  to  periodic  linear  motion 
is  illustrated  in  Figure  9.  Note  that  the  second  (outer)  link  of  the  pendulum  used  in  the  experiment 
was  slightly  shorter  than  the  first  (inner)  link,  hence  the  pendulum  is  slightly  different  from  the 
symmetric  pendulum  considered  in  many  of  the  previous  calculations. 

In  the  experiments  conducted,  motor  (and  consequently  oscillation  frequency)  control  was  open- 
loop.  Hence  we  only  make  qualitative  observations  on  the  stabilization  of  equilibria.  Using  this 
apparatus,  we  were  able  to  stabilize  all  equilibria  predicted  by  our  analysis  of  the  averaged  potential. 
These  equilibria  are  depicted  in  Figure  10.  In  conducting  these  experiments,  it  was  observed  that 
the  stabilization  sequence  agreed  with  the  numerical  calculations  of  the  last  section.  Specifically,  the 
hanging  folded  equilibrium  stabilized  at  roughly  the  same  motor  voltage  (^oscillation  frequency) 
as  the  inverted  folded  equilibrium,  and  the  fully  inverted  equilibrium  stabilized  at  a  much  higher 
motor  voltage.  Furthermore,  the  hanging  folded  equilibrium  was  observed  to  be  more  robust  to 
disturbances  than  the  inverted  folded  equilibrium,  and  much  more  robust  than  the  fully  inverted 
equilibrium. 

5.4  Basins  of  Attraction  for  the  Cart  and  Double  Pendulum 

To  this  point,  we  have  focused  on  the  stabilization  of  equilibria  and  ifurcations  of  equilibria  as  the 
result  of  variations  in  primarily  the  forcing  frequency.  What  arises  as  an  important  issue  for  the 
n-link  chain  that  is  not  so  important  for  the  single  link  chain,  and  is  important  for  the  type  of  open- 
loop  control  we  propose  in  general,  is  the  long  term  dynamics  of  n-DOF  systems  with  arbitrary 
initial  conditions.  The  issue  of  long-term  dependence  on  initial  conditions,  or  basins  of  attraction, 
is  approachable  for  SDOF  systems  because  of  the  existence  of  powerful  analytical  tools  such  as  Mel¬ 
nikov’s  method  and  invariant  manifold  theory,  and  essentially  because  it  is  possible  to  compute  and 
visualize  phase  space  separatrices.  This  analysis  has  been  performed  for  the  SDOF  rapidly  forced 
pendulum[18]  and  for  an  entire  class  of  SDOF  velocity  controlled  systems[15].  At  present,  the  theory 
for  higher  dimensional  systems  is  highly  specialized  and  much  less  developed. 

Our  long  term  goals,  for  which  we  provide  preliminary  results  here,  revolve  around  three  ques¬ 
tions: 

1.  To  what  extent  do  basins  of  attraction  of  the  averaged  system  describe  those  of  the  nonau- 
tonomous  system? 


Figure  10:  Stabilization  of  the  hanging  folded  equilibrium  (left),  the  inverted  folded  equilib- 
rium(center),  and  the  fully  inverted  equilibrium  (right).  In  conducting  these  experiments,  it  was 
observed  that  the  stabilization  sequence  agrees  with  that  predicted  in  the  previous  section. 
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Figure  11:  Basins  of  attraction  for  the  averaged  cart  and  double  pendulum. 


2.  How  can  forcing  or  other  physical  parameters  be  varied  to  favor  one  or  more  basins  of  attrac- 
tion? 

3.  How  can  the  averaged  potential  be  used  directly  to  estimate  either  the  basins  of  attraction  or 
the  probability  that  an  arbitrary  initial  condition  will  be  attracted  to  a  specific  equilibrium? 

In  addition,  a  detailed  understanding  of  the  basins  of  attraction  possibly  will  enable  to  quantify 
an  equilibrium’s  robustness  to  disturbances;  i.e.  quantify  how  much  of  a  disturbance  will  force  the 
system  into  the  basin  attraction  of  another  equilibrium. 

The  dimension  of  the  phase  space  for  even  n  =  2  makes  numerical  calculations  difficult  unless 
we  restrict  our  attention  to  some  two  dimensional  subset  of  initial  conditions.  One  reasonable  set 
of  initial  conditions  are  those  where  the  system  starts  at  rest.  To  develop  intuition  for  the  first 
question,  we  have  simulated  the  double  pendulum  system  with  small  damping  over  a  256  x  256 
grid  of  initial  conditions  which  cover  the  configuration  manifold  for  a  variety  of  inclinations  a, 
forcing  parameters  S  and  e,  and  inertial  parameters  fii  and  A^.  The  results  for  one  set  of  simula¬ 
tions  are  shown  in  Figures  11  and  12,  where  as  indicated  in  the  plot,  a  =  0.  In  the  figures,  we 
see  a  rough  correspondence  between  the  basins  of  attraction  of  the  averaged  and  nonautonomous 
systems,  which  is  typical  of  these  plots  for  other  parameter  combinations.  Noticeable  in  the  basins 
of  attraction  of  the  nonautonomous  systems  are  the  irregular  basin  boundaries  and  remnants  of  the 
regions  of  stochasticity  which  are  known  to  form  around  basins  of  regular  motion  in  periodically 
forced  Hamiltonian  systems.  The  basin  boundaries  represent  intersections  of  unstable  manifolds  of 
the  various  equilibria  with  the  initial  condition  space,  and  the  irregularities  suggest  that  at  the  value 
of  damping  used  in  the  simulation,  significant  tangling  of  stable  and  unstable  manifolds  persists. 
This  was  observed  in  the  case  of  the  SDOF  system  in  [18],  where  theorems  were  presented  which 
establish  bounds  on  the  tangling  of  phase  space  separatrices.  Also  note  that  the  folded  equilibria 
dominate  the  initial  condition  space,  as  suggested  by  the  results  of  the  previous  section. 
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Figure  12:  Basins  of  attraction  for  the  nonautonomous  cart  and  double  pendulum. 


6  Conclusion 

In  this  report,  we  have  considered  the  stabilization  of  a  cart  and  n-pendulum  system  on  an  inclined 
plane  by  high-frequency  periodic  forcing.  The  viewpoint  we  have  taken  is  that  we  wish  to  understand 
the  details  of  behavior  induced  by  a  very  simple  type  of  control  action  applied  to  a  reasonably 
complex  mechanical  system.  Specifically,  we  have  found  that  as  the  forcing  frequency  becomes  large, 
the  n-pendulum  system  aligns  itself  with  the  cart  motion.  These  results  were  obtained  by  deriving 
a  dimensionless  noncanonical  Hamiltonian,  and  then  averaging  the  noncanonical  Hamiltonian  to 
obtain  an  average  Hamiltonian.  The  results  followed  from  a  critical  point  analysis  of  the  averaged 
potential.  In  addition,  we  saw  that  the  results  obtained  through  the  use  of  the  averaged  potential 
may  be  connected  with  results  obtained  through  Floquet  theory.  Finally,  we  confirmed  the  results 
of  Section  2  for  the  vertically  forced  cart  and  double  through  a  numerical  study  of  the  system’s 
bifurcations  and  stability.  A  striking  feature  distinguishing  the  n~DOF  systems  from  the  more 
widely  studied  SDOF  case  is  the  sensitivity  of  the  results  to  nonuniformity  in  the  characteristic 
length  scales  of  the  system  components. 
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